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CHAPTER ONE 


INTRODUCTION 


This tutorial manual illustrates the programming features of 
APPLE (A Parallel Programming Language).* APPLE is based on 
K. E. Iverson's APL (A Programming Lanquage, Wiley, 1962), a 


language that uses generalized operators to concisely express 


mathematical algorithms on multi-dimensional data structures. 


The generalized operators of APL and its successor APIL\360 
permit a programmer to express manipulations of arrays almost as 
eaSily as he can express manipulations of scalars in 
conventional programming languages (e.g. FORTRAN, PL/I, etc.). 
The advantage of these operators is that when the programmer 
writes arithmetic expressions involving arrays, he does not have 
to go through the tedious process of writing nested loops to 
control the processing of the arrays. Instead, he is able to 
express the process as it conceptually peeves: in parallel on 


all of the elements of the arrays. 


That programmers tend to think in terms of parallel processes on 


i Ch Te < he <i <i ame Ca Paw a oe oro 


* APPLE 1s not to be confused with the RADC assembly language 
bearing the same name which was produced for an associative 
processor. 
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arrays is borne out by the APL\360 code produced by its large 
number of commercial users, who write their code as if it could 
be executed in parallel even though it is executed on a 
sequential computer. Many payroll evoavans: for example, treat 
the set of base eaiariee as a vector. Withholding taxes, 
voluntary deductions, and so forth, are then computed for the 
entire vector, at once, rather than on an employee-by-employee 
basis, because the same algorithm is used to orocess every 


employee. 


One of the primary benefits of APL has been the elimination of 
unnecessary loops and bookkeeping. For example, a programmer 
writes ioops far less frequently in APL in than in FORTRAN or 
PL/I. This tends to reduce the number of situations in which 


coding errors can be introduced into a program. 


However, because 4PzZ\360 1S not a complete programming language, 
it does not eliminate all such situations. The only control 
Operator in APL\360 is the GOTO operator. The basic control 
Operators--e.g., IF and DO--that have been included in FORTRAN 
and other languages dating from the 1950's ave absent from APL. 
Paradoxically, one can write elaborately eloquent arithmetic 
expressions in APL but must resort to the techniques of assembly 


language programming in order to perform them more than once. 


Dijkstra, Mills, Schorre and others have blamed a Ma JOESEY of 


programming errors on the unrestricted use of the GOTO 
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statement. Since one can control program flow only with GOTO's 


in APL, it was clear that there was a need for other control 
operators in APPLE if the possibility of coding errors was to be 


reduced Slgnificantly. 


Consequently, APPLE contains such features as conditional 
statements, conditional expressions, operators for writing 
loops, and case statements. These control operators eliminated 
any need for the GOTO operator. Therefore, there is none in 


APPLE. 


Other unique features of APPLE increase the clarity of 
exposition and simplify the coding process. These features 
further generalize APL's concepts and conventions for 


manipulating arrays and defining functions and operators. 


APPLE is sufficiently extensible that all of its operators can 
be defined in the language itself. A formal specification of 
APPLE is contained in the “ILLIAC IV Language Requirements 
Study: Final Report," SDC document TM-5074/000/00, 31 January 


1973. 
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CHAPTER TWO 


HOW TO USE THIS TUTORIAL 


In this tutorial, we assume that you have had some programming 
experience. We do not assume a sophisticated understanding of 
progranming languages. The manual is self-contained, so you 
should be able to learn how to use APPLE by reading the 
descriptions and working out the examples that have been 


provided. 


Because APPLE treats many mathematical concepts differently than 
do most other languages--certainly FORTRAN or PL/I--it is 
important that you read Chapters 3, 4, 6 and 7 closely, even 
though you may be familiar with many of the concepts. It is 
especially important chat you be aware of the differences in how 


to subscript arrays or evaluate arithmetic expressions. 


The tutorial is organized so that each chapter builds on its 


predecessors. 


Pay close attention to the discussions on the order of 
evaluation in Chapters 4 and 10. The concept is easily learned, 


but you must understand it thoroughly in order to program in 
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APPLE. 


APPLE has not yet been implemented on any computer. 
Consequently, we have not included any description of 
input/output or systems interface procedures. Those 
descriptions will be made available along with each APPLE 


implementation. 
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CHAPTER THREE 


WHAT IS AN ARRAY? 


The ILLIAC IV is a powerful computer. It was designed to 
data operands. For example, if we wanted to double each of a 


set of 50 numbers, the ILLIAC could double them all at once. 


In a conventional programming language, such as FORTRAN or PL/I, 
you would have to assign a unique name to each of the 50 numbers 
in order to write a program that would double each of them. One 
Simple way of assigning a unique name to each of the numbers is 
to declare an array that contains them. Then, each number would 
have a unique name consisting of the name of the array and a 


subscript (or index). The subscript would simply be a number in 


the range 1, 2,..4«,590. 


There are at least two ways to write a program that would double 
each of the 50 numbers. Ona way iS to Simply write 50 
assignment Statements, each of which sets a specific element to 
twice its previous value. Another way is to write a loop that 
will iterate 50 times and in which each element is replaced by 


its double. When coding a program, this latter alternative is 
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preterable since it requires less writing. Even more important, 
writing the loop reduces the possibility of your making 
keypunching errors. While neither of the alternatives in this 
approach seems very important, eonsider the problem of doubling 
100,000 numbers. In this case, we would have no alternative but 
to use an array and a loop since the program would be too long 


to write. 


While these techniques are perfectly acceptable ways of 
programming in FORTRAN or PL/I on a sequential computer like an 
IBM/370 or a PDP-10, we may ask: Why write a loop to double 50 
numbers on the ILLIAC when it only takes a couple of 
instructions in ILLIAC Machine Language? The answer is ehad 
these programming languages were not designed for computers like 
the ILLIAC, so there is no notation in the language to represent 
doubling all of the numbers at once since this is impossible on 


Ordinary computers. 


APPLE is not a conventional programming language. It is 
designed for use with computers on which you can double 50 
numbers all at once. In fact, APPLE is designed to run on a 
"computer" where you can double 100,000 numbers at once. (Since 
no such computer exists, the APPLE compiler makes the ILLIAC 
Simulate this imaginary computer. Thus, all you have to do to 
double the 50 numbers in APPLE, is to put them into an array and 


double the array. 
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Before we can show you how to write an APPLE program for 
doubling the array, we must first establish the terminology that 
is used to describe arrays and their properties. While all 
programmers are familiar with arrays, few programming languages 


treat them the same way. 
3.1 Vectors and Matrices 


A vector is a one-dimensional array of numbers. A vector is an 
ordered set of elements (i.e., a first element, second element, 
etc.), and we know how many elements there are. It does not 
matter how you write them down--horizontally, vertically, or | 
diagonally; the number of elements in the vector will not change 


and there will still be a first element, second element, and so 


on. The number of elements in a vector is called its 


eas, aioe ie eae <a aie dehy «aie alee Ely dae ee annie 


Example. If Y is a vector, then we represent its dimensionality 
by writing py. (9 is the greek letter "rho.") For 
example, if y is the vector consisting of the four 


numbers (14, 3, 2, 17), then pV is egual to 4. 


A matrix is a rectangular array of numbers. Each matrix has a 
nunber of rows and a number of cOlumns. Mathematicians can 
specify a particular element of a matrix by calling out, for 
example, the third element of the fourth row, or the fourth 


element of the third column. This is completely unambiguous. 
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They can also call out either the entire fifth row, or the sixth 


column, or one of the diagonals of the matrix. 


A matrix has two important dimensions: the number of rows and 
the number of columns it contains. The dimensionality of a 
matrix is defined as the vector whose first element is the 
number of rows and whose second element is the number of columns 
in the matrix. For a matrix M, the dimensionality of M is 
written pM. 
Example. If M is the matrix 
12 7 6 5 

8 10 22 15 


then oM is equal to the vector (2, 4). 


The dimensionality of a vector is defined to be the vector whose 


Only element is the number of elements in the vector. Thus, we 


‘Ge a aa ane Ge a dD ae aE CD Wty EE SRE emne «mney 


Suppose A is either a vector or a matrix. If the vector pA 
consists of only one element, then 4A iS a vector; if pA 

consists of two elements, then Ais a matrix. So we need a 
precise way of determining the number of Stenents in pA to 


decide whether 4 is a vector or a matrix. 


Since pA is always a vector, it makes sense to talk about the 
dimensionality of pA, i.e@., to talk about the one-element vector 
opA4 whose single element is just the number of elements in the 


vector pA. (Here, we write ppA to mean p(pA). The parentheses 
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are not necessary, sO we omit them.) It follows that if pp4 
equals 1, then pA contains one element, hence 4 is a vector. If 
eppA equals 2, then pA contains two elements, hence A is a two- 


Gimensional array (i.e., a matrix). 


It is cumbersome to talk about the dimension vector of the 

dimension vector of an array A. So we define the word rank to 
mean the value of the unique element of the dimension vector of 
the dimension vector of an array, 1.e., the value of ppA. Then, 


a vector is a rank-1 array, and a matrix is a rank-2 array. 


Mathematicians speak of row vectors and column vectors. In 
APPLE, these are not really vectors, but matrices. This is 
because a row vector always has two important dimensions: the 
number of elements it contains and the direction in which it is 
written. This.is also true of column vectors. In order to be 
consistent with the convention of listing the number of rows 
first, then the number of columns when we talk about the 
Gimensionality of a matrix. The dimension vector takes the 
following form. For a row vector Fk the first element of pF is 
always 1, the number of rows in the matrix, ana the second 
elenent is the number of elements in the row vector. Similarly, 
the first element in the dimension vector of a column vector is 
always the number of elements in the column vector, while the 


second element is always 1. 


It is possible to determine the number of elements in a matrix 
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by looking at its dimension vector. The number of elements is 
equal to the product of the number of rows in the matrix and the 
number of columns in the matrix. That is, the number of 
elements in a matrix is equal to the product of the two elements 


in its dimension vector. 


Example. If Mis a matrix having 5 rows and 7 columns, 
then pM equals the vector (5, 7). There are 
35 elements in ™ and 35 is the product of 5 


and 7. 


3.2 Rank-n Arrays 


In APPLE, the concept of an array is generalized to an arbitrary 
number of dimensions. We call this number the rank of the 


rray. For example, a rank~-3 array iS an arrangement of numbers 


if 


along the three coordinate axes of Euclidian 3-space. That is, 
the elements are arranged to form the lattice points of a 
rectangular parallelepiped. If A were such an array, then pA 
would be a vector (a, b, c), where a, b and c correspond 
respectively to the number of planes, poGs.and soins of A, and 
ppA equals 3. Similarly, there are rank-4 arrays, rank-5 
arrays, and so forth. It is easy to see that the number of 
elements in a rank-n array Ais the product of the elements of 


p As 


= 44 
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3.3 Empty Arrays 


APPLE permits you to work with an arrayA for which one or more 
of the elements of pAis zero. Since the product of the 
elements of pA equals zero, it follows that A contains no 


elements. 


You will occasionally have use for empty arrays. In fact, an 


empty array occurs in the following section. 


3.4 Scalars as Rank-0 Arrays 


A scalar is a number, as distinguished from a vector, matrix, 


quaternion, etc. A scalar corresponds to a geometric point. 


In APPLE, a scalar is an array that has no dimensions whatsoever 
associated with it. Consequently, there can be no elements in 
the dimension vector pS associated with the scalar5. This 
implies that pps equals 0, the number of elements in ihe wecke: 
oSe Since pps is the rank of the array S, we maintain 


consistency by calling a scalar a rank-0O array. 


a ED: 
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CHAPTER FOUR 


NOTATION FOR CONSTANTS AND VARIABLES 


In APPLE, a constant 1S a number the value of which never 


changes during the execution of a program. 


A variable is not a variable in the mathematical sense. Rather, 
a variable is the name by which you refer to a value that you 
wish to store someplace and access later. The value of a 
variable may change during the execution of a program, or it may 
remain constant. The Significant point is that the value of a 


variable can vary according to your needs, but the value of a 


constant is always the same. 
4.1 Notation for Numbers 


APPLE permits the use of integral and rational numeric 
quantities. These numbers are called integers and 


floating-point numbers, respectively. 


The precision of the ILLIAC permits the representation of 
integers n such that -2*8<n<2%8 (i.e. integers smaller in 


magnitude than 281,474,976,710,656). 


~ 13 - 
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Floating-point numbers are rational approximations to real 
numbers. The representation Panes for a floating-point number f 
is 2716384<¢,£4<216383, where the significant part of the 
mantissa is correct to 48 binary figures. The floating-point 
eéoreseneacien is automatically used for those integers that 


Cannot be represented in 48 bits. 


Integers are written the same way in APPLE as they are in normal 
mathematics, except that the negation sign is represented by a 
raised bar () so that it can be distinguished from the 
subtraction operator. Commas may not be used to separate three 
digit fields, because the comma iS an Operator that has a unique 


meaning in APPLE. 


Exanple. The number 1,234 is written 1234, while -50,762 is 


written 50762. 


Floating-point numbers are also written according to the normal 
arithmetic conventions. Here, too, the negation sign is used to 


represent negative numbers. 


Example. Pi may be written as 3.141592653583273, while 
244,397 as written 14.337 
You may also represent numbers in scientific notation, i.e., as 
the product of a number and some integral power of 10. The | 
‘mantissa does not need to be normalized. Here, the number is 


represented by writing the number, the letter #, and the 


integral power of 10 by which the number is to be multiplied. 


- 4 - 
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Exanple. The number -47335 can be written in scientific 
notation as either -4.7335x10*%, or 
~ -0047335x107 or-47335000000x107&. In APPLE, 
these would be, respectively, ~4.7335Fy0r 


~.,0087335EF 7OX 47335000000F 6 -« 
“.2 Notation for Variables 


Since one or more values is stored ina variable, we must have a 


means of referring to variables. We do this by giving the 


variable a name. 


A name consists of an alphabetic character followed by a 


(possibly empty) sequence of alphaneric characters. An 


coum amore ane Ge a ab ee a see <a Ce) ce Gite ite am Or aap ame 


a «me Wm see cm a de ep ae ED tp ae aS aie a a ae ab 


underscored digit, or an alphabetic character. A name may not 


contain any imbedded blanks. 


Example. The following are names: 


A 
Al 
Al 


4.3 The Equivalence Symbol 


APPLE uses the double-headed arrow (<”) to represent equality. 
This symbol is not an APPLE operator, but serves only as a meta-. 
linguistic device. Thus, when we wish to say that the content 


oe be 
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of the variable A is the number 3, we write A <> 3. 


If we want to be more precise and insist that 4 contains the 
scalar 3, as opposed to the vector (3), we would have to specify 
two facts: one related to the numeric value contained in A, the 
other related to the rank of A. In this case, we would write: 


A ee 3 
ppA <> 0 


If A had been the vector containing only the number 3, then we 


could have written either p4<«>1 or ppA +> 1, 
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CHAPTER FIVE 


BASIC ARITHMETIC OPERATORS 


APPLE provides the programmer with a large number of arithmetic 
Operators. These operators are designed to operate on arrays, 
rather than on scalars. Some of the operators; e.g., addition, 
subtraction, multiplication, division, exponentiation; are 
common to standard languages. The remaining operators are of 
the type commonly found in the mathematical subroutine libraries 


of major programming languages. 


The dperators are applied to entire arrays. The multiplication 
operator can be used, for example, to double all of the elements 
of an array without your having to write a loop. It can also be 
employed to multiply each element of one array by the 


corresponding element of another array. 


In this chapter, we will introduce you to each of the arithmetic 
operators and then explain how it works. We will subsequently 
describe how you form expressions involving more than one 
operator. Ina later chapter, we will show you how to 


generalize some of the arithmetic operators. 
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5.1 The Scalar Monadic Operators 


An operator is called monadic if it operates on only one 


argument (or operand). A scalar monadic operator is a monadic 


operator that is defined in terms of its effect on a scalar 


operand. 


Since each element of an array is a scalar, a scalar monadic 
Operator applied to an array operand A produces a resultant 


array B such that pA <> p& Each element of B equals the appli- 
cation of the operator to the corresponding element of A. 


A monadic operator is written to the left of its argument. 


In the remainder of this chapter, we will use the variables 


As By Ce: i, U, Vto represent the following arrays: 


oO Be OD 2 3 5 
A<>| 7 4 1 B+«+\7 4 1 
C <> (1.33, 1.33, 7.0, 0) D «> (2.72, 3.14, ~5.8, 148.3) 


U <> (1, 0, 1, 0) Y<«> (1, 1, 0, 0) 


where 


pA «> (2, 3) pB «> (2, 3) pC <> ( y) oD <> ( y) 
oU <> (4) oV <> (4) 


. ramet) 


The synbol + is used to represent the identity operator. For 


any array A, +A equals A. In symbols, we have +A +> A. 
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Example: 
rsp 2 os 0 
a 
5-1-2 The Negation Operator 
The symbol -1S used to represent the negation operator. (Note 


that - is different from the negation sign , which is only used 
for writing negative numbers.) For any array4, every element of 


ALS subtracted fromoO.,. 


Example: 


59.1.3 The Signum Operator 


The synbol ~ is used to represent the Signum operator. For a 
number x, Signum(x) is the function whose value is: 1 if x>0, -1 


if x<0, and 0 if x=0. 


Example: 


1 September 1973 2 8 System Development Corporation 
= 3 nf | TM-50747100700 


5.1.4 The Reciprocal Qpexrator 


The synbol : is used to represent the reciprocal operator. :B 
is defined for all nonzero arguments, its value is 1 divided by 


Be 


Example: 


3p oie eae Bere S Os 33030330333 
Oct DOS TIROS Ose-25 


hs © 
N> 


7 


The synbol * is used to represent the exponential operator. For 
any array A, *A equals e (2.7182818284...) raised to the 4 


power. That is, x4 is the natural antilogarithm of A. 


Bxample: 
eA <+/1,.353352832F 1 2.00855369221 1.000000000£0 
Gieercerrerer 1.831563889F 2 3.678794412F 1] 
5.1.6 The Natural Logarithm Qperator 


The synbol ® is used to represent the natural logarithm 
operator. For a strictly positive argument 8, ® is the 
logarithm of 3; to the base e. 

Example: 


@ 2 pa ana rey: 1.098612289 eee ereed 
4.9459101490 1.386294361 0.000000000 
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The synbol | is used to represent the floor operator. Fora 


number x, the floor of x is the algebraically greatest integer 


less than or equal to x. 


Example: If ¢C is the vector defined on page 18, then 


The synbol | is used to represent the ceiling operator. Fora 
number x, the ceiling of x is the algebraically least integer 


greater than or equal to x. 


Example: 


The symbol | is used to represent the absolute value operator. 


The absolute value of a number x is the algebraic maximum of x 


and -xX. 
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Example: 


5.1.10 The Random Integer Operator (Rol]) 
The synbol ? is used to represent the random integer operator, 
which is better known as "roll," as in the rolling of a die. 

The operand must be a positive integer array. For each scalar 
n, the result is a normally random integer selected from the set 
{1,2,eee,n} Or {0,1,.2-¢,n-1} according as the index origin is 7 
Or 0, respectively. (See Section 6.1 for a discussion of Index 


Origin.) 


Example: 


5.1.11 The Logical Complementation Operator (NOT) 


it once eee Te eeelioammeal ‘esto canta NN weettiee eens enmee coe wee) Cae eRe mans semeet wnt einen enmmee andiiD ew nits stant eenaenra ee exe me que eR wee 


The synbol ~ is used to represent the logical complementation 
Operator. The operator is defined Only on the set {0,1} and 


transforms 1 into 0 and O into 1. 


Example: 


1 September 1973 System Development Corporation 
TM-5074/100/00 


The symbol !is used to represent the generalized factorial 
operator. ! is not defined for negative integers. For all non- 
negative integers n, the result of applying this operator is n! 
If x is not an integer, the result of applying this operator is 
the gamma function applied to x+1. If n is a negative integer, 
n! is undefined. 

Example: 


'C <> (1.188192811, 4.08546585, 5040, 1) 


The synbol 9° is used to represent the operator that multiplies 


its operand by 


Examples: 


OC <> (4.178318229, 4.178318229, 21.99114858, 0) 
5.1.14 Summary of Scalar Monadic Qperators 


The 13 scalar monadic operators are summarized in Table I. 
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TABLE I 


SCALAR MONADIC OPERATORS 


Symbol Name Meaning 
+ Identity +A <> A 
~ Negation -A «> 0-A 
41a21£f A>0 
x Signum xA | O if A=0 
“4 if A<o 
+ Reciprocal tA <> 14A 
k Exponential *A «> eA 
@ Natural @A<+> In 
Logarithm 
[ Floor Algebraically greatest integer < B 
[ Ceiling Algebraically least integer 2 B 
| Absolute |4 <> { A if A20 
Value -A if Ao 
? Random Inte~ Random Integer between J@QRG and A 
ger (Roll) 
~ Logical Com- ~A +> 1-A (for Ae{0, 3} 
plementation 
t Generalized 14 <> \ A factorial if A <> |LA 
Factorial (TA+1) if A#LA 
O n Times OA <> 1XxA 


5.2 The Scalar Dyadic Operators 


An operator is called dyadic if it operates on two operands. A 


terms of its effect on a pair of scalar operands. 


If Aand B are two arrays such that pA +* p2, then the elements 


of Aand Bmay be paired according to their positions in the two 


arrays. We say that an element from Aand an element from® are 


CUD A CUED NED 6 aie ee GE 


within their respective arrays; that is, the subscript that 


identifies the one element also identifies the other. The 
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application of a scalar dyadic operator to two such arrays 4 and 
B produces a resultant array © where e¢ «> pA +> PB and an 
element of © corresponds to the result of applying the operator 


to the corresponding elements of A and ®é. 


A dyadic operator is written between its arguments. The variables 
A,B,C,D,U,V used in the discussion of dyadic operators 


are defined on page 18. 


9-2.1 The Addition Operator 


rh Sa ag GD ner ait GERD CD aie al eae a Ch a ae a ae een 


The symbol + is used to represent the addition operator. 


Example: 


9.2.2 The Subtraction Qperator 


a ae aw CD MD ae eh EE Cane cane ee ae em eae ee Th Ga Ce ae om en 


The symbol - is used to represent the subtraction operator. The 
argument on the right of the operator is subtracted from the 


argument on its left. 


Examples: 


‘ 


A-B o> (-u 0 5 
0 78 7 
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The symbol xis used to represent the multiplication operator. 


Example: 


> am om ie > aa a dae a ath ee 7 aun 28> Ge «a a ate 


The synbol + is used to represent the division operator. The 
argument on the left of the operator is divided by the argument 
on its right. The operation is defined for nonzero divisors 


only. 


Example: 


The synbol | is used to represent the residue operator. If m#0 
and n are numbers, there exists an integer q such that n = mg + 
r, where 0 $ r < {m{.~ The symbol r represents the residue of n 
modulo m. The definition that follows is extended to cover the 


case m = 0: the residue of any nonnegative n, modulo 0 is equal 


to n, but remains undefined for n < 0. 


1 September 1973 System Development Corporation 
TM~5074/100/700 


When we write 4A/|/i, we mean the residue of 2 modulo 4. 


Example: 


A|B <«>70 O 5 
O 0 O 


CLD ee (0.06, 0240, 14.2, 2853) 


dan a aT en CLD akg ce ete «ne Cita -—_ CR GM > Gis hte Ge 


The synbol | is used to select the algebraic minimum of its two 


operands. 
Example: 
ALB <*+[°>2 3. 40 
ea, 
CLD «>» (1.33, 1.33, 5.8, 0) 
9e2-/ The Maximum Operator 


The synbol [ is used to select the maximum of its two operands. 


21° = 
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Example: 
Alb «7272 3 °5 
\7 44 
ChD «> (2.72, 3.14, 7, 148.3) 
5.2.8 The Exponentiation Operator 


The symbol * 1s used to represent the exponentiation operator. 


In order to raise 4A to the B power, you write Ax38. 


Example: 


AxB <> UL,OD00000000LZ0 2.,.700000000F1 0.00000000080 
8.23543000075 2.,.560000000F2 1.000900000H0 


R*xeA <> 


H0000000r 1 22 7FO00O000ORT 1.00000000070° 


2%0 
8.235430000F5 3.906250000#F 3 1.000000000F0 


DeC +> (3.784222315, 0.2183149959,  2279798.4168, 1) 


Note that 4*p is not always defined. For example, A*0.5 is the 
square root of 4, which is defined only for nonnegativeA. Of 


course, —32*0.2 +» "2 since 0.2 <* t:S5and §2*5 «<> ~32, 
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5.2.9 The Logarithm Qperator 


The synbol ® is used to represent the logarithm operator. The 


logarithm of , to the base B is written BeA. By definition, A 
and B must be strictly positive and we may have A +> 1 if and 


only if B «+ 1 The common logarithm of A is written 10@A, 


Examples: If 


T <>/10 10 10 
10 10 10 


T@B <+/(0.3010299957 0.4771212547 0.6989700043 
0.8450980400 0.6020599913 0.0000000000 


5.2.10 Lhe Circular Function Operator 


The symbol © is used to represent the family of operators for 
all of the trigonometric and hyperbolic functions, which are 
collectively referred to as the circular functions. This 
applies to the trigonometric functions since they are defined in 
terms of the unit circle, and to the hyperbolic functions as a 
consequence of the relations Saale = i sin z, cosh iz = 


i cos z, and tanh iz = i tan z, where 12=-1. 


A circular function is invoked by writing AOB, where the value 


of 4 is used to identify the particular circular function, as 


99 = 
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follows. 


The following trigonometric functions are defined for angles in 
radian measures: 
40B is equivalent to sing. 
20B is equivalent to cossB. 
30B is equivalent to tang. 
“4090p is equal to arcsin Bp, where 12|B 
“20B is equivalent to arccos B, where 12|B 
“30B is equivalent to arctan 2B. 
Three functions are useful in trigonometric identities: 
uoBp produces the principal square root of 1+£2, 
0OB produces the principal square root of 1-£2, 
where i12|B _ 
~4OB produces the principal square root of 1+B2, 
where i1<|2B. 
(Here, we abused the language somewhat by writing 52 to mean 
Bx2e This was done to avoid confronting you with a complicated 
equivalence like ~40B «+ ( 1+B*2)*0.5 at this early stage. See 
Section 5.5 for a discussion of the priority of operators and 


APPLE's bracketing conventions.) 


The hyperbolic functions: 
50B is equivalent to sinh 6 
60B is equivalent to cosh 8B 
70B is equivalent to tanh 8 
~50B is equivalent to arcsinh 82 


“60B is equivalent to arccosh 8, where 221 
-“70B is equivalent to arctanh B, where 1>18, 


5.2.11 The Logical Conjunction Operator (AND) 


The symbol 4A is used to represent the logical conjunction 


operator, AND. This operator is defined Only on the set {0,1} 


wt “BO ass 
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and is completely defined by its action on the vectors Uand / 


(see Section 5.1): 


UKRY ee C1. Ory. 04. 0) 


The synbol vis used to represent the logical disjunction 
Operator, OR. This operator is defined only on the set {0, 1} 


and iS completely defined by its action on the vectors JU and V/V : 


5.2.13 Tbe Exclusive Disjunction Operator (Not Equal) 


The synbol * is used to represent the exclusive disjunction 
Operator. This operator's domain is extended to the set of real 
numbers; its range is {0, 1}. A#B +> 1 if and only if A and B 


are unegual. 


Example: 


A#B <«>{1 0 1 
cee ae 


UzV <> (0,1,1,0) 


a. 24 = 
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5.2.14 The Equality Qperator 


The synbol = is used to represent the equality operator. The 
domain of this operator is the set of real numbers, while its 


range is {0,1}. 4=2 «+ 1 if and only if A and & are equal. 


Example: 


U=V <> (1,0,0) 


ie wile a ab aakene eu ~—_ ai ilies ee 18S se tm 


The synbol * is used to represent the NAND operator. NAND is 
defined to be the logical complement of AND. The domain and 
range of NAND are {0, 1}. NAND is defined by its action on the 


vectors J and /: 
UN: > Cg Te Ls. 1) 


5-22.16 The NOR Operator 


ae a ae ae i feet ae ah ae a aa wa 


The symbol »¥ is used to represent the NOR operator. NOR is 


defined to be the logical complement of OR. The domain and 
range of NOR is {0, 1}. NOR is defined by its action on the 
vectors {/ and /: 


Uw~V «> (0, 0, 0, 1) 
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The synbol < is used to represent the less-than operator. It 


maps the reals onto {0, 1}, A<B+> 1 if and only if A is less 


than i. 


Example: 


5.2.18 [he Less-Than-or-Egual Qperator 


The synbol * is used to represent the less-than-or-equal 


operator. It maps the reals onto {0,1}, A<8 «> 11f and only if 


A 1s not greater than 86. 


Example: 


5.2.19 The Greater-Than-or-Equal Qperator 


The synbol 2 is used to represent the greater-than-or-equal 
Operator. It maps the reals onto {0, 1}. A2B +> 1 if and only 


if ais not less than 2. 


= 93. 4 
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Example: 


The synbol >is used to represent the greater-than operator. It 


maps the reals onto {0, 1}, 4>8 «+> 1 if and only if A is greater 


than A 


Example: 


Jimi Heh: OU 
0 0 9 


5.2.21 The Generalized Combination Operator 


ma <> a a ee ab ae ie ie i ae Ch SP Gre tie Gin ap Gh TN tee abe ae ie aie ES CREB wee CE cues 


The synbol - is used to represent the generalized combination 
Operator. 4! is the number of combinations of 2 objects taken 4 
at a time. If a and b are nonnegative integers, the number of 
combinations of b objects taken a at a time is given by 
C(b,a)=b!/a!(b-a)! The generalized combination operator uses 
almost the same formula, but replaces the factorial operator 
with the monadic generalized factorial operator. Consequently, 
the generalized combination Gpeeaeae is defined for all 
arguments for which the generalized factorial operator is 


defined. 
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Example: Let 
he Oy: By My, “Date Ge. 0) 


Then 


5.2.22 The Random-Selection-Without 


The synbol ? is used to represent the random-selection-without- 
replacement operator. The result of writing A’7 (where A<B) is 
a vector », such that pk «+ A and the elements of # are randomly 
selected without replacement from the set {1, 2,..., Db} or {0, 
1,e--, b-1} according as the index origin is 1 or 0, 
respectively. (See Section 6.1 for a discussion on index 
Origin.) The operator ? is defined only for nonnegative scalar 
integer arguments. The operator can be used to simulate the 


dealing of bridge hands, for example. 


Examples: 
678 «> (8,7,3,2,4,5) 
678 «> (8,5,2,3,1,4) 


£3762 eo (O54 74104 2 oad 54 by TO 5 12s 94 oe LO) 


ae te 
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5.2.23 Summary of Scalar Dyadic Qperators 


The 36 Scalar dyadic operators are summarized in Table II (see 


next page). 


« 36: 
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TABLE II 
SCALAR DYADIC OPERATORS 
Symbol Name Meaning 
+ Addition A+B 
‘ Subtraction A-B 
x Multiplication AxB 
: Division A+B 
| Residue A|B «> B(mod_ A) 
[ Minimum ALB «> mint A, B} 
[ Maximum AIB «> max{A, B} 
ze Exponentiation 4*«B «+ 48 
® Logarithm A@®B +> log |B 
Seen 0 see Circular Domain 
—703 arctanh B 1>|B 
“603 arccosh 8B B21 
~50B arcsinh B 
~4uOR (~14+B*2)*,5 1<|B 
“302 arctan B 
“208 arccos B 12/B 
~40B arcsin B 12|B 
OOB (1-Bx2)*.5 12|B 
10B Sin B 
20B COS B 
308 tan B 
HOB (1+B*x2)*.5 
50B sinh. B 
60B cosh 8B 
70B tanh B 
A AND AAB 
<j OR AVB 
2 Exclusive OR A#B «+ (AVB)AAw®B i1f A, B ec {0,1} 
Inequality A#B 
= Equality A=B 
w NAND A®B <> ~AAB 
» NOR A¥B <> ~AVB 
< Less Than A<B 
< Less Than Or A<B 
Equal 
> Greater Than A2>B 
Or Equal 
~ Greater Than A>B 
' Generalized A!B «+> (!B)#(!4)x!B-A 
Combination 
"4 Deal 
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5.3 Right and Left Identities 


Suppose that the symbol * is some scalar dyadic operator. If 
there is a number L such thatl*B <> B for every value of 8, L 
is called a left identity of x. Similarly, if there is a number 


R such that A®R «> A for every value of A, then Ff is called a 


identity £, it follows from elementary algebra that R <+ £Z. 


Table III summarizes the identity elements of the 36 scalar 


Gyadic operators. 


a BG te 
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TABLE III | 
IDENTITY ELEMENTS OF SCALAR DYADIC OPERATORS 


Operator Left Identity Right Identity 
+ 0 0 
- none ¢) 
x 1 1 
+ none 1 
x none 1 
@ none none 
| 0 none 
O none none 
Vv 0 

A 1 1 
"~ none none 
~~ none none 
t 1 none 
[ —c = co 
L 00 co 
? LORG LORG 
The following identity elements apply only to the domain {0,1}: 
> none 0 
> none 1 
< 0 none 
< 1 none 
= 1 1 
a 0 0 


5-4 Right-Associative Operators 


An operator xis associative if, for any 4A, 2. G we always have 
Ax(BxC) +> (A®B)xCG For such an operator, there is never any 
ambiguity; you can write AxBxC and everybody knows what you 


Meane 


However, not all operators are associative. For example, the 
subtraction operator is not associative since, e.g., (5-4)-3=-2, 


while 5-(4-3)=4. (The first interpretation is called a left 


= 39 =, 
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associatign; the second is a right association.) It is ambiguous 


to write 5-4-3, since it is reasonable to interpret this 


expression with either grouping of terms. 


Before we tell you which interpretation APPLE makes, let us look 
at a slightly more involved expression involving subtraction, 


say a-b-c-d-e-f-g. 


First, let us look at the parsing (((((a-b) ~c) -d) -e)-£)-g. 


Since -b=(-1)b, we have 


(((( (a-b) ~c) -d) ~e) -£)-g = ((( (Cat (-1) b) + (- 1) c) + (-1) dd) + 
+(- 1) e) + (-1) £) + (-1) 9 
at+(-1) (b+c+tdtetftg) 


a> (b+c+dtet f+gq) 


Since addition is associative. We see that this parsing is 


equivalent to subtracting the sum of all the other terms from 


CP? 2S Gn Gm cme in czar GD Gi Gab ae GS Gate aE eo a 


parsing. 


If we had used the other parsing, we would have 


(a-b) +c- (d- (e- (£-g) )) 
(a~b) + (c-d) +e- (f£-g) 
(a-b) + (c-d) + (e-f£) tg 
(atctetg) —- (b+d+f) 


a- (b- (c~(d- (e~ (£-g))))) 


That is, you take the sum of the first, third, ... terms and 
subtract the sum of the second, fourth, ... terms. This is 
called a right=associative parsing. In a sense, the right- 
associative parsing of a-b-c-d-e-f-g is more interesting than 


the left-associative parsing, a-(btctdtetft+g). Division and 


exponentiation are also nonassociative operations. They are 
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interpretation than when given the left-associative 
interpretation. For example, with left-association Ax«B«CxDxFE is 
just Ax(BxCxDxE). Under right-association, it is equal to the 


more familiar 4x(Bx(C*(D*E))) from algebra. 


Since right-associative parses are generally more interesting 
for nonassociative operators than left-associative parses, all 


ee ae aye 


APPLE Operators are treated as right-associative operators 


This applies not only to expressions involving repetitions of 
the same operator, but to expressions involving mixtures of 


APPLE operators as you will see in the next section. 


5.5 Bracketing Conventions and Operator Priorities 


Consider the expression 2+3x4. The rules of algebra say that 
this expression evaluates to 44, 1.e., to 2412. This is because 
algebra assigns a higher priority to the multiplication operator 
than to the addition operator. If you had Wank this 
expression to evaluate to 20, you would have had to parenthesize 
the quantity you wanted evaluated first, writing (2+3)xu 
instead. (Of course, you would have removed any possible 
ambiguity by writing 2+(3x4) when you wanted the expression to 
evaluate to i4, but this is not necessary when you know the 


operator priorities.) 


With a few exceptions, most programming languages follow the 
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standard operator priorities of algebra. However, there are 
some differences in the way some Languages treat an expression 
like i12:u4x3. In some languages, the result is 12+12, while the 


expression evaluates to 3x3 in others. 


So far we have covered 49 APPLE Operators. It would be 
difficult for anyone to remember the relative priorities between 
such a large number of operators. Many of the assigned 
priorities would appear artificial. There are approximately one 
hundred operators in APPLE, hence the problem is non-trivial. 
So in order to simplify the problems of learning APPLE, there 
are no Operator priorities whatever in the language. Instead, 
the right-associative parsing convention is extended to 
expressions involving a mixture of Operators. If you want an 


operator to take priority over some other operator, all you have 


to do is parenthesize that operator and its operands. 


Returning to the expression 2+3x4, we see that it is equivalent 
to 2+(3x4) , that is, 14. But 4x3+2 is dauiipatene to 4x(3+2) or 
20. Hence, expressions are not necessarily commutative in 
APPLE. If you wanted to have 4x3+2 equal to 14, you would have 


to write either (4x3)+2 Or 24+3x4. 


Consider 12+4x3. The right-associativity of + and x means that 
the expression ig equivalent to 12:(4x3), ie-e€.1. The 
expression means that 12 is to be divided by whatever is on the 


right of the + operator. The quantity on the right of + is 4x3, 


ee UD a 
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i.e. 4 multiplied by whatever is on the right of <x. Since that 


is just 3, we see that we are dividing 12 by4x3, i.e. by 12. 


Note that monadic and dyadic operators can be intermixed in an 
expression. For example, 4+-30:5+6>7xeg is a complicated- 
looking expression. Let us add parentheses according to the 
right-associativity convention: 4+(-(30:(5+(6>(7x(@8))))))- (We 
know that the minus sign is an operator since "minus thirty" 
would have been written ~30. You cannot interchange § and - 
since 30 has an eeeeee only on the magnitude of "thirty," while 
-3+5+6>7x8 changes the sign of everything to its right. We know 
that - and @ are monadic operators since each is preceded by 


some other operator rather than an operand.) 


Now to evaluate the expression. We start with the most nested 
subexpression. eg is the natural logarithm of 3, i.e., 
2.07944154178nNd7x2.0794415417 +> 14.5560907919. Next, 
6>14.5560907919 ++ 0. 5+0 +> 5, and 30:5 «+ 6- Next, “6 «+ -6 


and u+ 76 <> ~9e SO we see that 4+-304%5+6>7x@8 «> 2. 


Expressions are never evaluated backwards in APPLE, although 
they are evaluated from the right. When you write 4-5-6, the 


result is 5 Since 5-6 <> “1 and 4-1 +> 5. 
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CHAPTER SIX 


ELEMENTARY ARRAY MANIPULATIONS 


In this chapter, you will be introduced to a class of operators 
that are useful in manipulating arrays. These manipulations 
include the familiar process of extracting One Or more elements 
from an array by subscripting. Since those elements that have 
been extracted are arrays, recall that scalars are rank-0 
arrays, you will be forming a subarray of the original array 
each time you subscript into it. The subarray may consist of 


more than one element. 


Subscripting is only one means of forming subarrays of an 
Original array. You will be introduced to techniques for 
forming subarrays consisting of elements satisfying some set of 
properties, as well as techniques for taking various cross 


sections of an array. 


We will also describe ways of rotating and transposing arrays, 
of combining several arrays to make a bigger array, and of 


changing the dimensionality of arrays. 


es A 
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6.1 Index Origin 


IN) 


ubscripting is the process by which you specify one or more 
elements of an array. The subscript of a specific array element 


is known as the index of that element. 


There is some confusion between the way programming languages 
refer to the first element in a vector. In some languages, that 
element has an index of 1, while its index is 0 in other 
languages. Depending on the programmer's particular needs, one 
of these indexing conventions is often preferable over the 


other. 


Since only you, the programmer, know which indexing convention 
is preferable for your personal application, APPLE leaves the 
choice up to you. The index origin is the value of the index of 
the first elements of a vector. The index origin is contained 
in the rank-0 array JORG. iI0ORG normally contains the value 1. 
If you want to specify its value, you begin your program with 


either 


or 


The value of JORG will remain constant throughout the body of 
your program. For the time being, the values of JORG are re- 
stricted to 0 or 1. Eventually, the language may be extended 


a eee eels enn 


to permit arbitrary integral values of JORG. 
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6.2 Indexing of Arrays 


In the following subsections, you will be given the necessary 


vocabulary and notations for subscripting arrays. 


tebe Onn amet CRD at <> 


Recall that if X¥ is some array and pX <> 0, then X is an empty 
vector (since xy contains no elements, and ppX¥ «> 1, SOX isa 


vector) . 


We will be using the empty vector often enough to require giving 


it a name for easy reference. The name of the empty vector is 


Suppose you wanted to generate a vector that consists of all of 
the permissible values of indices, in ascending order, for some 
vector VY. Clearly, such a vector consists of pV elements. If 
IORG «+ 0, then this vector is (0,1,..., ite’) if JORG +> 1, 


then this vector is (1,2,...,pV). 


In APPLE, you need only write ip V to produce this vector. The 
Operator 1 produces the desired vector. In fact, all you need 


to do to get a vector of length n, where n20, is write in. 


= 6 = 


1 September 1973 System Development Corporation 
TM-5074/100700 
Example: If JORG «> 1, then 


15 «<> (1,2,3,4,5) 
11 <> (1) 


while if 7oRG <> 0 then 


Lo. es (O54552, 364) 
10 ++ (0) 


If we write 0, from the definition of 1, we should get a vector 
of length 0. A vector of length 0O can only be the empty vector 
€~ Hence, regardless of whether JQ&&G «> 1 or [ORG +> 0, we will 


always have 10 «> e€. 


So far, iN is defined only for nonnegative integers WN, whereN 
is either a scalar or a one-element ‘vector. We will soon extend 
the definition of 1 to cover all rank-0 and rank-1 arrays 


consisting of nonnegative integers. 


~ Ey = 
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If you want to indicate the k-th element of some vector VV, where 
1SkSpy, the appropriate index would be k+lORG-1, For example, 
if you want the first element and JORG «> 1, you want the index 
to be 1+1-1, i.e. 13; while if JORG «> 0, then you want the index 
to be 1+0-1, i.e. 0. Correspondingly, if you wanted the fifth 
element, the index would be either 5 or 4. You should convince 


yourself that kt 7ORG-1is always an element of ipVvV. 


Scalar subscripting is exactly like subscripting in other 
progranming languages. In order to select the k-th element of 
Ve you write V(K), where K is the appropriate index 
corresponding to JORG and k. For example, suppose J@fG +> 1 and 


V «> (1,.5,7.9,3,4,1)* Then v[i] <> 14 VI3] «+ 7, VE61 <> & If 


When you subscript a vector with a scalar, the result is a 


scalar. 


You are not restricted to using scalars as subscripts, however. 
If you subscript a vector with a vector, the result is a vector 
of the same dimensionality as the subscript vector. For 
example, if A <> [2,1,5,4] and Vis the vector we used in the 
previous example, if JORG +> 1, VLA] +> (5,1,3,9) and JORG +> 0, 


then [A] <> (7,5,4,3) 
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When you subscript a vector with an array, the cesult has the 


same dimensionality as the subscripting array. 


Example: 
If 7oRG ++ 1 and j is the rank-2 array 


= 
4 
+ 
IN 
me OD 
b> O71 
NOmes 


then 


ViW] +> 27 3 


6.2.4 Subscripting of Arrays 


The nunber of elements in an array A is the product of the 
elements of the vector »4. Since A has ppA coordinates, any 
subscript of A must be composed of ppA components. These are 
separated from one another by the delimiter (;). The first 
subscript you list applies to the first coordinate of A; the 
second one applies to the second coordinate, and so forth. If 

K «+» LORG+k-1 and JS the k-th coordinate of A, it is 

reguired that J be in the range JORG< I < IORG + (pA)EK]-1. 
That is, a coordinate subscript must lie in the range 1,...,(pA)Li] 
if) fORC = OF Oransy COALKIe? 22 20OnC <s Oo. Wiis 

is always equivalent to saying that a coordinate subscript is an 


element of the vector ipA[K]. 


eo: = 
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Example: | 
If forng «> 1 and Ais the rank-2 array 


5 
8 7 2 
Then pA <> ( 
The subscript for the first coordinate must be either 1 or 2; 
the subscript for the second coordinate must be 1, 2, 3 or 4. 
When you write A[i;i], you specify the scalar 3. We also have 


AL 228) #* 2,AL 2.4) «> © and AL2s2) e785 


Just as you can subscript a vector with an array, you can also 
subscript an array with an array. The result is an array 8% 
whose rank op is equal to the sum of the ranks of the 
coordinate subscripts. The dimensionality of 2 ,pB, is the 
vector that is composed of the dimension vector for the first 
coordinate subscript, followed by the dimension vector for the 


second coordinate subscript, and so forth. 


For example, we can subscript 4 with the vector (1, 2) for the 
first coordinate and the scalar 3 for the second coordinate, 
thereby producing the array composed of A[1;3] and A[2;3]. The 
result must be a rank-1 array Since pp(1,2) «> 1 and pp3 <> O 
and 1 + 0 +» 1. The dimension vector of the result is the 
vector (2, €) «+ (2) Since p(1, 2) «> 2and p3 «> e€ . Therefore, 


we must have 4A[1 2; 3] «> (5, 2) 


APPLE permits you to write vectors in subscript expressions 
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either with or without parentheses. Of course, for the sake of 


Clarity, you could have written 4[{(1, 2);3] if you had wished. 


Now, suppose we have a matrix 


so that pC «> (2, 3). When we write4A[c;3 4 ] , what should the 


result, », be? 


We know that ppp «+ 3(since popB «+ 2 and pp(3,4) «<> 1), and 

pi «> (2,3,2) the catenation ( composition) of pc and p(3,4)). 
For Simplicity, call i} «> (3,4). We can deduce the elements of 
2, aS follows. P[1;1.1] must correspond to the’.1:1] -th row of 4 


and the | 1] -th column of 4A. Hence, 


Dit ie). => iis ee 5 
Similarly, we obtain: 
Pi t2tseo) es ALL eae) <>. 6 
DLi3;23:1J <7 AL1:3) <> 5 
DiL1ises2) <> ALIsy] =>.-6 
Pi o3e4 ) ee hls ad ee 2 
YL13;33;2)] <% AL2;4] «> 1 
PL2s13;i1] <* AL2;:3] «> 2 
PES eo } = 2224) eer 7 
PL2:2:1] <> AL1;3] «<> 5 
DL2:2:2] «+ AL13;4] «> 6 
DL 22321) => AL2s3)] <> 2 
Di223320) <e- ALI a) => 4 
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We can graphically represent D as 
5 6 


6 
1 


ho & 


D«+«--> 

2 1 

>. 6 

a 
Finally, suppose you want the first plane of 2). You could write 
Di13;1 2 33; 1 2). The result would be a rank~-2 array of 
dimensionality (3,2} (Why?) It would consist of D[1i;1;11], 
2 Slee eae and so forth, as expected. Recalling that 13 +? (1,2,3) 


and 12 «> (1,2), you could also write 2f13;13312], 


By convention, instead of writing 1(p4)[K] asa 
subscript for the k-th coordinate of an array 4, you 
can elide* the subscript for that coordinate, writing 


any required semicolon separators as you normally 


would. 


Formally, we have the equivalences 


AlsdsK3.0e3¥iZ] «<> ALi pAVLLORGI373k pee ae 
ALS tek tee so + ALL: lo )LAtLORG J 3K 3% et 2 
A Pele hee ctl AL IidsK:...3¥310pA)CIORG+ 1+ppAl] 


Hence, instead of writing Dl1;133;12], you can write 211334, and 
the APPLE compiler will deduce the content of the elided 


coordinate subscripts. 


* —Elide--to omit. 
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6.3 The Ravel Qperator 

when we constructed the array 2? in Section 6.2.4, the order in 
which we listed its elements was significant. We started with 
the element whose subscript was 7oRG in each coordinate. 
Subsequently, we allowed the right-most coordinate subscript to 
vary most rapidly, then the subscripts in the coordinate field 
second from the right, and so on. This ordering is called an 
would follow had they been placed on the individual wheels of an 
automobile odometer. (Wheel xk of the odometer is numbered with 


the elements of i(pD)[K], starting with JORG.) 


It is sometimes useful to view an array aS a vector. In APPLE, 
the comma (,)is used monadically to represent the ravel 
operator. The ravel of an array is the vector whose elements 
are those of the original array in the odometer order. In 


particular, the ravel of a scalar is the vector whose only 


element is the scalar. 


Thus, we see that 


opi 49° 0556.5 5629.1 424155504241) 


6.4 Array Index Generation 


Here, we generalize upon the definition of the monadic operator 


1° Suppose 4 is some array and WV <> pA, so that W is a vector. 
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Then, we define ,4 to be the matrix of dimensionality — 


p(iN) <> 


((0,4),(ppA))= That is iN has asS many rows as there 


are elements in 4 and as many columns as the rank of A. Thus, 


ywcontains a column for every component of a subscript of A and a 


row for every element of A. 


The rows 


is the index associated with the element 


to the i-th element of ,A (i.e., 


(,A)(I]). 


Example: 


of iW are in odometer order, so 


Let us refer back to the array 
Section 6.2.4. oD <> (2,3,2) 


We have 
1 i 
i ee 
1 2 
1 2 
13 
| 4,3 
wD <F 94 
21 
2. 2 
2 2 
2. 2 
2.3 


NRO FR RN RD eR DD FD Fe 


(iW)CIIJ 


 f 


that the i-th row of iW 
of A that corresponds 


is the index in A of 


D constructed in 
so let us look at ipD. 


Now let us look at a few examples of the concordance 


between ioD and D. 


First, 


NO 
> 


D<> 
2 
<) 
| 2 
and ,D «> (5,6,5 


4 
6 
A 
26, 
and (1pD)(3;] <> (1 


2,1,2,1,5,6,2,1). Now,(.,D)€3] <> 5 
2,1)- We see that p[1;2;1] +> 5- 


> 
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We also have, (;,D)(11] «> 2, (19D)L113;] «+ (2,3,1) 
and D[2;3;1) +> 2. 


In this example, we assumed thatror¢g «> 1. If 
IORG <> 0, it follows that 


[200000 
Dm rFRFROO 
RBOrFORO 


1pD <> 


bs ps Pe PS Pd pS 
NMEE OO 
FPOnNrRPHO 


6.5 The Subscript Generator 


Suppose Ais an array and JI is some row of 190A. There is an 
element of 7; for each component of 4. But in order to subscript 
A by 7, there must be semicolon delimiters present between the 


elements of J. 


The monadic operator ;/ is used to generate scalar subscripts. 
Its only action is placing a semicolon between the elements of a 


vector. 


Thus, if pA <> (4,7,3.,5,8), the vector IT «> (2,3,2,1,6) is 
certainly a row of the matrix ipdA. We have 


AL3/Z£] «> AL23;33231;6]- 
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6.6 Partial Subscripting 


Suppose Ais an array and K is an element of ippA, that is, K 
is an element of (0,1,..., 1+994) 1£ IORG «+ if OF an element of 
(0,1,..., 1pp4) if JORG «+ 0. Suppose that I is some array in 
which each element of I is an element of i(pA)([K] so that Iisa 
valid subscript array for component K of A. Then, if the only 
component of A that you want to subscript is component kK, you 


may do so by writing AL[K]I]. 


K must be an integer scalar, or the integer scalar content of a 


variable. 


Example: Suppose JORG <> 1 and A is the rank-3 array. 


FMR 
OnwW Wm 


OF EO 
OOF Wo 
On WhD 


Then pA «> (2,4,3). Now, [1;;] is the first plane 
of A. This could be written [[1]1], which 

says the subscript 1 is to be applied to 
coordinate [1] of A only. | 


When we write A[L2]2 3], 
have pAL[2]2 3] <> (2,2,3) 


2: 4.3 

| 07 6 
ACC2]2 3] <> 

4.23 

4 0 7 


We can also look at A[[3]2 3]. 
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Peres, pl Tati -as > Ce dhe op oanea 


Al C3J]2 3] «> 


6.7 Reshaping Arrays 


Any array can be transformed into a vector by use of the ravel 
Operator. Any array , can also be transformed into an array 8B 
of different dimensionality. To do this you use the reshape 
operator op. The reshape operator is dyadic, while the shape 
Operator p that yields the dimension vector of an array--its 


shape--is a monadic operator. 


When you want to transform 4 into 8 so that pB +> #, where # 1s 


some vector of nonnegative integers, you write B <« Pod 


8 must contain a number jy of elements equal to the product of 
the elements of fF Since »B ++ RF, B receives its elements from A 


according to the formula: 
BL3;/CapB)CII] «+ AL3/(CipA)LZLORG+i |I-LORGI) 


This formula is another way of saying that a vector y is 
constructed from enough repetitions and partial repetitions of 
the elements of (|, that VY contains as many elements as BP will 


contain. In odometer order, the first element of B will be the 


G7’ x 


1 September 1973 System Development Corporation 
| TM-5074/100700 
first element of JV, the second element of /! will be the second 


element of VY, and so forth. 


Example: Suppose JORG «> 1 and Ais the rank-3 array 


& 2 3. 4 
5 6 7 8 
9 10 11 12 
A «> 
dd: 2k 5-26 
Ly 18 LS 20 
21. 22 23 24 
Then pA +> (2,3,4)6 
Recall that 
on. ee Oe ee ea Oe De boc tat Oy io Oe Oe eae hee eee Oe 
If i24 were reshaped so that its dimension were(2,3,4), 
it would be equivalent to A. That is, 


(2?,3,4)p124 <> A, 


The vector i60 has 60 elements. But, when we write 
(2,3,4%)p160, we are only uSing the first 2x3x4 «> 24 
elements of i160. Therefore, 


(2,3,4%)9160 <> (2,3,4)p124 +> A. 


Verify that 


NO O71 
wo 
~J 
= 


C25 3s4) 017 <> 


Gd 
Pg 
Cn 
oD) 


Note also that 
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1 2 3 #4 
5 6 7 8 
@ 10> ia. ta 


(6,4)pA «> 


4 5 6 
- 8. SS? “20 2k ole2 


eS 
NO 
Ww 


(4,6)pA «> 
12 14 15 1617 18 
19 20 21 22 23 24 


Finally, (10)o93 is the scalar 3. This is true since 


0(10)3 +> §& and pp(10)93 «> pe «<*> 0 
6.8 Arithmetic Array Manipulations 


This section introduces the very useful classes of operators 
that enable you to perform arithmetic processes on the elements 


of an array. 


6.8.1 Vector Reduction 

It is often desirable to obtain the sum or product of the 
elements of a vector. In APPLE, you write +/V in order to 
obtain the sum of the elements of V ; you write x/V in order to 


obtain the product of the elements of/V . 


The sun of the elements of a vector is usually considered to be 
a scalar. In APPLE, while V is a vector, +/V’ is a scalar. 
Thus, the rank of +/V is equal to the rank of Y/Y reduced by 1. 


For this reason, +/V is read as the "plus reduction of /." 
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Example: 
ue ings 6 
K/ (O23 Ste 5 
You can use any of the operators from Table III (see Section 
5.3) in conjunction with the slash to form a reduction operator. 
The effect of reduction is to place the dyadic scalar operator 


between the elements of the vector operand, and then to evaluate 


the resulting expression. Consequently, -/(2,3,4,7, 5,2)#*2-3-4-7-" 


“414, If pV «> 1 for any operator #, */V +> (10)pY. 
If pV +> 0, for any Operator ®, %/V <> (10)o7, where J is the 


identity element associated with ee 


It is useful to note that x/pA 1s the number of elements in the 
array Ae When A is a scalar,pA «> € .« Since x/p4 is always 
equal to the number of elements in any array, regardless of its 


rank. 


Note that we always have {|/i «+ lorRG for any le and 
[/iN «> N-14+I0RG-e- We will use this notation frequently when we 


<n 


write subscript expressions that are independent of the value of 


LORGs 


6.8.2 Array Reduction 


Suppose A is a rank-n array and suppose *® is a scalar dyadic 
operator from Table III (see Section 5.3). Then */[K]A is the 


application of x over the elements of coordinate K of A. The 
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result is an array of rank -1. Its dimension vector is 


obtained from pA by suppressing (9A)([K]. 


Example: Suppose JORG <> 1 and pA «+> (2,4) where 
> 7 10 
4 3°41 2 
Then 
+/LijJA <> 
+72 lA. => 


( 
( 
-/Li]JA <> ( 
-/L2]A +> ( 


Bo+<> 1:3. 28: oS 


61 = 
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Then 


+/{1] 


+/[2] 


+/{3] 


aes: 


Peee 


7 Eos 


Leas 


L/C2] 


Eva 


B 


B 


B 


B 


B 


b 


> 


4--} 


<> 


<-> 


+--> 


<--> 


30 
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When X is the last coordinate of A, you may elide the [[/1ippA4] 


and write *®/A. 


Hence, 


if A is the matrix from the example, we 


roe as 
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have 


+/[L2]A <> +/A 


If A is any one-element array (i.e@., x/pA <> 1 , then #x/[K]A is 
a one-element array of rank Of i+tppA. For example 


+/(1,1,1,1,1)93 «> (1,1,1,1)93 while+/(10)p3 <> (10)p3. 


The accumulation operator \ is the analogue of the reduction 
operator. If x is one of the operators from Table III (see 
Section 5.3), #\V is the x -~reduction of the elements of /, 
Starting with the i-th element of vy. In other words,x\y isa 


vector of the "partial sums" obtained in the evaluation of x/y. 


Example: Suppose 
V «+> (5,4,3,2,1) 
then 
x\V «> (120, 24, 6, 2, 1) 
and 


As in the case of array reduction, *\[K]JA is the application of 


over the elements of coordinate x of A and opx\[K]A es AG 


Example; If JORG «> 1 and M is the rank-2 array 


x 63:< 


x 
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M<--> 9 1 10 2 
Then 

+\C1]M <*> 9 1 #10 2 
13 11 10 3 

+\L2]M <> 4 13 12 2 
g 1 7 3 

F\L1JMéM <> 9 4 10 2 
7 7 7 3 


[\€2]M@ «+ 9 2 "2 °2 
6.9 The Catenation of Vectors 


If V and Ware any two vectors, then you can construct a new 
vector x¥ «<> V,Ww where pV,W <> (pV)t+pW and x[ipV] <> V 


Xi CpW)o (pV)+i1pW] <> W 


That is, the first py elements of y,yware the elements of V and 


the last pW elements of V,W are the elements of i. 


V,W is called the catenation of y andy. The catenation 
Operator(,) is a dyadic operator. If either y or yw is a scalar, 
it is treated as a vector, so that y,w is a vector of dimension 


(p,V)+o.We 


Example: If Uis the scalar 1 and Vand W are the vectors 


V «+> (2,3,4) 
Woe (5,657,548 ) 


Then 
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VW Ah LOS Nee hb 7 ek) 
U,W «> (€1,5,6,7,28) 
U,V,W «+> (1,2,3,4,5,6,7,8) 


6.10 The Interval Operator (Optional on First Reading) 


The synbol / is used to represent the monadic operator that 
produces an interval vector. The argument of Jis a vector 
LEN,ORG,S.¢d LEN,ORG,S is the vector of length Lil’? whose least 
element equals ORG. S must be either 0 ori. If S++ 0, 
successive elements increase by 1; if otherwise, the elements 


decrease by 1, 


Note that J is independent of the value of ORC. 


6.11 The Subarray Function (Optional on First Reading) 


It is often desirable to work with a subarray 2 of an array A 
where ppB «> ppA and the hyperplanes of 8B are adjacent 
hyperplanes of A. For example, if pA <«* (3,5,7) we might want 
to construct an array 


B +> A[T1 234 3 233 4 5] «+ Ald 2,1,03s/ 3,2,13¢ 3,3,3] 
The dyadic operator 4 is useful for this purpose. The right 


argument of Ais 4 and the left argument is a rank-2 array F 
where pF «> ((pp4),3)-e The rows of F are the vectors required 
by the interval operator / in the expression above. [hat is, 


the elements of # satisfy the formula 
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Bo«> FAA<*ACHY PLIORG: 13d PUCLORG+131;...3¢ FEE /ipp4A3]) 


Thus, for the example, F would be the matrix 


e tT. 0 
ee Se 
3 og: OS 


6.11.1 The Whole Array Qperator (Optional) 
The noOnadic use of A On an array A produces the matrix F such 


that FAA <> Ae F satisfies the following conditions: 


of oF (99A),3 
FO3L/13] +«*% (pA) 
FL31+b/13] +> (ppA)eLORG 
FL3P/13] +* (ppA)po 


Example: 


If pA «+ (2,3,5,7) and JORG +> 1 


Then 
21 0 
AA <*> 3 1 Q 
5: “i: © 
qd. 2 8 


If opB +> (3,5,6,8,2) and LORG +> 0 
Then 


Ab «<-> 


RO Ow MH OM WwW 
COO CO © 
Oocd’coo 
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6.11.2 The Cross Section Operator (Optional) 


An array cross section is obtained when all of the component 


subscripts are either scalars or elided. For example, if A is a 


rank-4 array, the following are some possible cross sections of A. 


A 

AL1i3s33] 

AL e323 J 

ri es eee 
Ate 3 34] 


The dyadic cross section operator A is primarily used for 
formalizing the subscripting of an array by scalars. The right 
argument of A is the array to be subscripted. The left argument 
of 4 is a rank-2 array G, pG «+ (ppA),2. The elements of the 


first column of G are either O or 1 as follows: 


Tf coordinate K is to be elided, then G[LK;] «> (0,0) 
Tf coordinate K is to be subscripted by the scalar °, 
then 


GLAs) eo C128) 


Example: If pA +> (2,3,5,7), IORG +> 1, and 


amme 


OF ORM 


Then 


GAA <> AL23;343;6] 
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6.12 Compression and Expansion 


aan oe oe we a cD ein a i aie eae ae 


that pU +> pX, we can form the subvector JU/xX (read "the J 
compression of xX"), where pU/X <*> +/U. The elements of U/<X are 


the ALIJ such that VII] «> 1, 


Example: If U and V are the vectors 


then : 
U/V «> (1,5,11,13) 
Since 
~U «<> (0,1,1,0,1,0,0) 
CMU) 7 Vee 025 347) 


6.12.2 Phe Logical Compression of an Array 


The logical compression operator is extended to arrays, as 


follows. Let A be an array, and let U be a logical vector such 
that for some component J of A, pU +> (pA)[I]. Then, U/[I] A 


the U compression along coordinate I of A is defined as 


UIGTA A ee ALLE? UY eon 


a Be pe 
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Example: Suppose J@0RG +> 1 and Aand JU are 


i] <7 C14140:) 


1 
uy 
7 


Co UT RO 
m 


10 11 12 
A «> 13 14 15 
16 17 18 


is ee © mee 
22 23 24 
25 26 27 


Then 


u/t1] A + 


and 


U/C2] A << 13 14 


and 


U/€3] A «> 10 11 12 


6 9- 
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\ 


If 7 <» [/ippA then U/[I] A may be written as U//, 
6.12.3 The Logical Expansion of a Vector 


The logical expansion of a vector is the analogue of the logical 
compression of the vector. In this case, if VY is any vector and 
U is a logical vector, where pV <> +/U, then J/V is the vector 

having zeros wherever JU has zeros, and whose remaining elements 


are taken in order from V/V. 


Example: 


(4,450 51460) VCl a 2g 8) oO 01425058 0) 


6.12.4 Ehe Logical Expansion of an Array 


The logical expansion operator is extended to arrays, as 
follows. Let 4 be an array and let JU be a logical vector such 
that for some component J of 4A, +/U <> (pA)LIIJ . Then U\LIJIA, 
the ’ expansion along coordinate J of A, is defined as the array 
whose dimension is given by 


eCUNEL IAI LL S670 
and for Xan element of 1094, p(U\LIJA)CE] «> (pA)LK] 


where XI and for every element / of 190J, 


CU\NA)LCU TIS] «> ACL II]+/ULidg J xl) 


= FO: = 
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The following example clarifies the situation. 
Examples: Suppose JOFG +> 1 and Uand Aare 


U ce C120525:0)) 


2 4 
6 8 
A <-> 


1. 33 


Then 


U\LIJA «<> 


and 


U\L2]A <> 


U\E3]A aa og 


If 7 «+ [/ippAe then y\[I]A may be written as U\A. 
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6.12.5 The Relation Between Expansion and Compression 


The following relation occurs between the expansion and 
compression of any array A by any logical vector JU. 
U/LIJU\LIIJA +7 A 


6.13 Prefix and Suffix Vectors 


A prefix vector is a logical vector P whose first # components 
are ones and whose remaining (pP)-R components are zeros. A 


Se ssanadlidinmsiesteaainnsadiaemessetiiten icsgilincaenedl CR ir i aes et 


ones and whose remaining (o0S)-# components are zeros. 


The dyadic operator for specifying a prefix vector is 4% the one 
for specifying suffix vectors is w In both cases, the left 
argument Lis the length of the vector to be constructed, and 
the right argument Ff is the number of ones required in the 


vector. 


Example; 


Prefix vectors are useful for selecting the first * components 
along some coordinate of an array; suffix vectors are useful for 
selecting the last#? components along some coordinate of an 


array. 


=» Jo < 


1 September 1973 System Development Sl parr 
TM-50747100700 


Example: Suppose /0RG +> ltandA is the array 


1 2 3 4 5 
6 7 8 9g 10 
11 $42 13 14 195 
16 17 18 19 20 


Then 
| 1 
( lig ? )/L1 JA Sate 5 7 Q Q 


C4w2)/LIIJA <* 45 46 17 18 


(su3)/L2]A 


4 

v 
— 
Gd 
> 
ty. é 
jd 
C3 


6.14 The Monadic Transpose Operator 


If 4 is any rank<-n array, where n22, the monadic transpose of 4, 
written of, is equivalent to 4 with lits last two coordinates 


interchanged. If A is a matrix, then &4 is the transpose of 4. 


Formally, for (ppA)22, 
pA «> pALri 2tppA), (f/ippA), 141 /i1ppA4) 
and for any row L of ip%A, we have 
AL;/Z] «> (&84)03/L0M]] 
where 
M +> (1 2+ppA), (T/ippA), 1+1/1p04) 
Example: If 


1 2 3 4 
5 6 7. 8 


9 10 11 12 
13 14 15 16 
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Then 


QA «> 


6.15 The Take Operator 


The take operator is a generalization of the prefix and suffix 
Operators. It permits you to form a subarray pg of an array 4 by 
writing 7tA where 7is a vector of integers such that p7 <> pod 


and (|7) s$ pAe. 


The effect of the take operator is that 


pB «+ |pA 


and coordinate J of B receives the first 7[7] elements from 
coordinate J of A if 7[I]20; otherwise g8 receives the last 7/7] 


elements of coordinate J of A when 7[I]<0o. 


The take operator can be expressed, as follows. 
tA <-> FRA 


where 


Fo«> Q(3,ppA)PCIT), CCC pT) pLORG)+(7<0)x(pA)-17),(pT)po 


APPLE has a convention that makes it easier to write certain 


take vectors. If you Only want to subset the last few 
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coordinates of A, you may write 7tA where (p7)<ppA. This is 


equivalent to writing (((~(ppA)ap7)/pA),7)tA. 


Also, if you want to subset a few coordinates of A and keep the 
rest intact, you may write 7t(IJA where J is the vector of 
coordinates in which you are interested. This is equivalent to 
writing tA, where p22 «> ppA and BLI] «+ 7, while 


(l41pp4)/B «> (L#1pp4)/pA.~ 


Examples: Suppose J¢RG +> 1 and ppd +> 3, where 


1 2 3 4 
> 6 FF « 
og). Ae it. 2 


Then 

(i,2, 3) tA +> 5 6 7 

(1,2, 3)4A «> 6 7. 8 
(1, 2, 2)+A4 «> 23 24 

6 7 8 

_ 40. 3 49 

(2, 3)tA +> 18 19 20 

D9. D3. Ou 


Since, by convention, (2, 3)tA «> (2,2, 3)tA 


Finally, 


- 7S. x 
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an 


| 6 7 8 
(1, 2)t£1,2JA «> o 40 44. 49 


5 6 7 8 
9 40. 44. 49 
~24(2]A «> 
17 18 19 20 
24.29 93. oH 


6.16 The Drop Operator 


The drop operator , is the analogue of the take Operator. When 
you write +A, the first or last |7[I] components of coordinate J 


of A are suppressed according as 7[I] 18 positive or negative. 


Fornally, 
TYA «> GAA 


where 
G «> Q(3,pp4)p((pA)-17,°0CC (CoP) plLORG)+007),(ppA)po 


The conventions mentioned for the take operator also apply to 


the drop operator. 


Examples: Suppose JORG «> 1 and A is the array defined in 
Section 6.15. 
Then. 


Cie. 3)d4 > Clg is 1je3 


2¥A +> 
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6.17 The Reversal Operator 


If Ais any array, then $[7]A is the reversal of coordinate 


I of A. Formally, 


oCIJA «> HAA 


where 


H +> Q(3,ppA)p(AA)DL3ZORGI,(AA)L31+L0RGI,((ppA)prI)=1ppA 


If I «+ [/ippA, then I may be elided. 


Example: As in Section 6.15, suppose JORG +> 1 and 


A <> (2,3,4)p124 


Then 
io 2H 15 2.6 9 10 11 12 
17 18 19 20 5 6 7 8 
21. 22.23.24 41 2 3. 4 
OL1jJA +> O(2]4 +> 
1 2 3. 44 21 22 23 24 
5 6 7 8 17 48 19 20 
9 10 11 12 13 14 15 16 
4 3 #2 1 
8 7 6 5 
i211 10 9 
oA <> 
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6.18 The Mask and Mesh Operators 


The mask and mesh operators combine the power of the compression 
and expansion operators in an extremely useful way in that they 
permit you to construct a new array from the logically selected 


portions of the other arrays. 


Let U be a logical vector and suppose X and Yare arrays where 
oX «> p¥. Then, the y mask of coordinate x of x and y is 
written U/(KI](X;Y). The resultant array is composed of the 
elements of coordinate K of Y, corresponding to the ones of J, 
and the elements of coordinate K of X, corresponding to the 


zeros of JU, 


Formally, 
opU/CKI(X3;Y) <> pX «> oY 


U/CK] U/CKICX3sY) «> U/CKIY 
(~U)/CKI] U/EKICX3Y) «> (~U)/CKIX 


Example: Suppose JQRG «> 1 and 


1 3 5 7 2 y 6 8 
A <> 9 1113 15 B<- 10 12 14 16 
17 19 21 23 18 20 22 24 
Then 
2 3 5 8 


(1,0,0,1)/CL2](A3;B) «+ 101113 16 
, | 18 19 21 24 
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6.18.2 The Mesh Operator 


Let U be a logical vector and suppose * and / are arrays where 
opX «> ppYe Suppose that for some component J, we have 

(CCppX)pl)#1ppxX)/px +> (CopY)#ippY)/py. That is, except for 
coordinate [, the dimension vectors for xX and Y are equal. [In 


addition, 


(oX)CL] «> +/~U 
(oY) LI] «> +/U 


Then we define U\(7IJ(X;Y), the U mesh of coordinate I of X and Y 
as the array of dimension 

pU\ETICX3¥) *> ((I-LORG)+(pX)),(pU), (141-L0RG) 4 (0X) 
That is, the resultant array consists of the elements of the 
U expansion of coordinate I of Y and the ~U’ expansion of 


coordinate J of xX. 


Example: Suppose JORG «+ 1 and A and B are 


4° > 3. Cy 100 200 300 400 
A<* 5 6 7 8 Bo ++ soo 600 700 800 
910 11 12 


1 2 3 Ly 

Then 5 6 7 8 
(0,0,1,0,1)\LiJ(A3B) «<*% 100 200 300 400 

Qo “40: 12. 12 

500 600 700 800 


6.19 The Rotate Operator 


Suppose 4 is an array and Wis a scalar integer. You can rotate 


a O° x 
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the elements of coordinate I of A cyclically N positions to the 
left if wWeo, or ii positions to the right if WV is negative. This 
is done by using the dyadic operator 6, writingW¢A . A cyclic, 
left rotation means that the left-most elements migrate around 
to the right-most position in their coordinate. Since the 
coordinate is of length (pA)[I] the rotation moves elements 


(oA)LIJ|N positions to the left or right. 
Formally, NO[IJA «> ALLIIJ(C(CppA )pLORG)+(pA)L III CC pp4 pW) +i (pAt(T) 


Example: If fORG +> 1 and 


Then 


7O,2jA «<> 19[2]A+> 16 17 18 


22 23 24 


and 
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2oC1JA<> 16[1]A +> 4 


If Ais an array such that ppA <> 1+p99B8 and 
(C(ppB)pL)4ipp£)/pl «> pA(i-e., the dimension vectors of A and P 
are in agreement after (pB)[I] is Suppressed, Ao[lr]B 1s defined 


as the rotation of coordinate 7 of B by 4. 


Formally, if 1 is a row of ipB, then 


(AOLTIB)L3/L] «> (AL3/C(ppB)pl)#1ppB)/LIOFAB)LC] 


where ¢C <> (((poB)pl)=1p0B)/L 


Fo <> Q((((pp8)pl)=19pB)/(031)),°00 (p08 )pl)=a1ppB)/(03L) 
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Example: Suppose JORG <> 1 and 


RO NO 
ho 
G2 RO 


A <> 


Then 
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AOC SIR ae 10-42 15 


and 


AOlLiIB <> bh 5 25 


AolijB «<> 15 13 14 


6.20 The Catenation of Arrays 


In Section 6.9, we defined the catenation of vectors. Two 
arrays A and g may be catenated along coordinate J provided 


either: 
(1) ppA +> pps 


or (2) 1 «<> |(ppA)-ppB 


a Oo: S 
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or (3) 8B or Ais a scalar 


and (4) the coordinates along which A and B are to be 
joined are of the same dimension. 


The meaning of condition (4) is the following. 


(a) If (ppA) «+ (pp&), then 
((ppA)plI)#1ppA)/pA +> (CC ppB)pl)#ippB)/pB 
That is, pA and pf are identical for every 
coordinate, except possibly coordinate [. 


(b) If 1 +> (ppA)-ppb, then 8B is considered as if its 
dimension vector is (((ppA)pI)=1ppA)\(1392). This 
dimension vector must satisfy condition (a). 

(c) If 1 <> (ppB)-ppA, then A is considered as if its 
dimension vector is (((oppB8)pl)=1ppB)\(13;pA). This 
dimension vector must satisfy condition (a). 

(d) If Aor Bb is a scalar, it is treated as though its 
dimension vector is identical to that of the 


nonscalar in all components, except for a unit 
component. 


Then the catenation along coordinate J of B and A, written 
A,lijJ8, is the array cC such that 


CCOT11( pA) CII] ]«> A 
CCLIIJ(CpADLTI+(pB)CIU <> £ 


Example: If JORG +> 1 and 
2 4 6 8 
Beri0 12 14 16 


Then | 


i: QU os 
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And a1 ees 


ON WOH 
c | 

SCs cain 
ee) 


1 3 5 7 2 4&4 G) 8 
A,C2]B <-> o:-44. 29 45 “40 “22. 14 2G 


If 
50. Si 2-55 
Bu 55 56. 57 
C <> 
58 59 60 61 
62 63 64 65 
Then 


A,£l1]¢ «> 50 51 52 53 


po 
Go 
an 
~] 


A,l1jJ2 «> 


tO 
ay 
ped 
fe 
Ww 
= 
qn 


7 i @ Ss <7 2 
A,l2]2 <> a bd Ae Aes 2 


6.21 The Lamination of Arrays 


When two arrays, 4and 4, are laminated together on coordinate 
I, a new coordinate that has indices i2is formed before 
coordinate 7. The argument4 fills the first index of the new 


coordinate; the argument 4 fills the second index of the new 
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coordinate. The notation is 4:(I2J2 where either pA +* pBorA 
ox , iS a scalar and [I is an element of i11+((pp4)[ppB)I£ A or B 
is a scalar, it is considered as if it were reshaped to have the 


sane shape as the nonscalar argument. 


Example: Suppose IORG «> 1 and A and B are the same as in 


Section 6.20. 
Then 


ALAS so 


ASL 218 => 


10 12 14 16 


A:L3]B <> 


A:[1]5 «> 


On 
nn 
on 
on 
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6.22 The Dyadic Transposition of Arrays 


The dyadic transposition of an array is a generalization of a 
monadic transposition of an array. The VY transpose of 4 is 
written Vi, where Vis a vector containing one element for each 


coordinate of A. 


The values of the elements of V indicate the dimension of the 

resulting array. If you want the result # to be of rank pe V 
must contain, in any order, at least one of each of the elements 
from ipph Further, the elements of Y are limited to the values 


contained in ipphk 


If the same value appears nore: than once in’, which must occur 
if (poR)<pV then the repeated dimension of fF is to be formed 
from more than one dimension of 4. For example, suppose that 
pepA «> 4, JORG «> 1,and V «> (2,1,2,2). The result & <> V&A 
must be an rank-2 array since 2 «> [/V. The first dimension of /?? 
is formed from the second dimension of A. The second dimension 
of fF is formed from the first, third and fourth dimensions of 4. 
Thus, the elements to be selected from A are of the form 
All;J;I;I], where 7 and J are scalar integers. This is the 
diagonal passing through A[1i;1;1;1] and through the first, 
second and fourth dimensions of A. The diagonal contains no 


more elements than the shortest of the dimensions from which it 


ls taken. 
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In summary, the vector V of positive integers must satisfy the 


conditions: 


qy) oV +> ppAe There must be an element in V for 
every component of A. 7 


(2) Every element of VY must be an element of ippA. 


(3) Every element of if /V must be an element of’. 
(For example, if ppA «> 6 and the largest 
element in yis 5, V must contain all of the 
elements of i5. Since VY must have six 
components, one of the elements of 15 must be 
repeated.) | oo 7 


Then, the transpose of 4A by V is defined as: 
(a) ppV&A <> 1401 /V)- 


(b) For each element J Of 1pp V&A , 
o( VALI] <% L/CV=(pV)pT)/pA 
(c) For each row £ of ipV%A, 
(VQA)L3/L «<*> AL3/ECVIIA 
In order to better understand dyadic transposition, study the 


Following detailed examples. 


Examples: Suppose JORG +> 1, 

pA «> (5,7,3,8 Je 
and y «+ (2,1,2,2). 
Then if f > VQAe 
ODL <> De 
From property (b), we see that 

(oF) Ca ete 1/002 942 291 st) 705375058) = 170051,0,0) 7053759858) SS 7 

(pR)L2I] <> L/(C2, 13235 2 = (1, ding 1% ys CS. Tay 0). S def Cogen tl) Se 2 


Therefore, (ok) «> (7,3) 


From property (c), we can determine the > 
mapping between elements of Aand RA. 
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(VQA)E131] <> ADL3s/(1,1)02 1 2 23) <> AL13;13;131) 
(VRA)L132) «> AL3/(1,2)02 1 2 21) <«% AL231325;2] 
(VQA)£C13;3) «> AD3/(€1,3)02 1 2 23) «> AL33133;33] 
(VQA)C23;1] <> AL3/(2,1)02 1 2 23) «7% AC1323131)] 


(VQ@A)E53;2)] +> AL3/05,2)02 1 2 21) «> AL2353232] 

(VQA)C733] +> ALCS/(7,3)02 1 2 2]] «> AL3373333) 
Thus, you see that the elements of V84 +> AlIs/;i;I], 
where J 1S an element of i17and J is an element of 13. 
Suppose we want X «> W8A where W <> (3,1,2,2) Then pp <> 3 © 
and oX ++ (7,3,5)- The mapping between 4 and yYis given 
by 

XCIL373KJ) «> ALK3I3d3d] 


where 7 is an element of 17, ¢ is an element of 13, and 
K is an element of 15. 


Finally, suppose 


Verify that the following are true: 


(see next page) 


Pa: ae 
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(1,1,1)0B «> (1,13) 
(1,1,2)8B «+ 13 14 15 
1 4&7 10 


(1,2,1)8B <> 13 16 19 22 


1 
(2,1,1)0B <>» 5 17 
| 9 


C25 252) 80 ee 13 27 21 


C2yie2 OR + 


on £ pe 
= 
a 


1 16 
CBe25 1.00 =e. 2 a7 
3.18 


(1,2,3)8B <> B. 


1 4 7 10 
2 5 8 11- 
3 6 9 12 
(14332788 <> 
13 16 19 22 
14 17 20 23 
145 18 21 24 
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nies: 


GO Ne 
> 
a xg 


uo16 
17 


(35152208 a Sad 


1. 23 
4 16 
4 19 
LO 222 


2 14 

7 
(3,2,1)QB «> 8 20 
1 


23 


oO 5 
6 18 
SAE 
42 24 
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Yoe5 66 
16 17 18 
(24142 )9R: <> 


10 11 12 
22 23 24 


1 44 7 10 
Lo Ae EO: a2 


(2.7.4 )N8 «2 «8 8 4d 
14 17 20 23 


3 6 9 12 
15 18 21 24 


~~ 9) = 
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CHAPTER SEVEN 


EXPRESSIONS, STATEMENTS AND PROGRAMS 


So far, you have seen a number of APPLE's arithmetic and 
manipulative operators. You have nearly enough information to 
write a program. However, you need to learn about APPLE 


Statements and expressions before you can write a program. 


In this chapter, we will cover everything you need to write a 
Simple program uSing APPLE operators. In the following 
chapters, you will be introduced to operations on arrays, as 


well as how to define your own functions and operators. 


7.1 Elementary Definitions 


a te tate LT ES 


operands. Expressions always have a value. 


A statement is an expression whose value is the empty vector. 
Statements always perform some action. In FORTRAN, they 
include, for example, assignment statements, DO statements, IF 


statements, and so forth. 


mG F 
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7.2 Contormability Conventions for Scalar Dyadic Operators 


Let us reexamine the simple expressions 4A®B where* is some 
scalar dyadic operator. Until now, we have stated that this 
expression is defined only when pA +> pB. Two such arrays are 


said to be conformable. The result is an array a where 


oC <> oA «+ pB and for any subscript Z from ipC, we have 


Clh3/L] <> AL3/LI*BL3/f] 


We will now extend the definition of conformability with respect 


to a scalar dyadic operator. 


First, suppose only one of the two operands, operand 4A, 1S a 
one~element Sere: (Aisa one-element array if and only if 

1 «+ x/ode Thus, 4 might be a scalar or a one-element vector.) 
We subsequently define A* to be the array C +> ((pB)pAxB , (If 
B was the one-element array, we would have C +> A¥(pA)p8 That 
is, the scalar operand is applied to every element of the non- 


scalar operand. 


Examples If JORG +> 1, then 


1+15 ¢%(2,3,4,5,6) 


If 
“4 3 y 
A «> 2°79 8 
Then 
2 6 8 
Ax 2 «> “4 14 16 
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If both operands are one-element arrays, then AxB is defined to 
be the array c where C +> ((ppA)I ppb)p(,4)#,% That is, the 
rank of the resultant one-element array equals the maximum of 


the ranks of the operands. 


If neither of the operands is a one-element array, the arrays 
are conformable only when the two arrays satisfy one of the 
following two conditions. 
(1) Assume that (ppA)>pe8 and (-ppB)tpA +> OB That is, 
the last ppF elements of pdAare identical to oB. 
(If (ppB)>ppA,everything works when you mentally 
interchange eheiee-nsness) Then, the arrays are 


conformable and C +> Ax®(pA)pB 


(2) If (1) is not true, but there are non-negative 
scalar integers Mand W such thatM<ppA and V<ppoBb 
and 


M+opA «> NYO 


Then A and B are conformable only if MtpA +> IMp1i or 
N+pA «> Hpi. The result is of rank ppC <> (ppA)looB, 
pC <> (O(ppC)p(opA),De1)I OC ppC)p(op8),2Pp1 anda 
C +> (CpC)pA)#(pC)pB, where D +> |(ppA)-ppB. 
In all other cases, A and PB are non-conformable arrays. 
Examples: If pA +> (2,3,5) and (98) «> (1,1,1,2,3,5) then A and B 
are conformable and pAx#®B +> (1,1,1,2,3,5) If 
pf «+> (1,4,3,2,8) and oF <> (1,1,6,4,3,2,8), then £ and F 


are conformable and e#H#F <> (1,1,6.4,3,2,8). If 
OG + 0147534540428) and pf eo Cot Pd 24 28s 
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Then G and ii are conformable and 
oGHIL «> C424 24 ie 24 Og eee Go) Ss 


If oV «> (1,2,3,4) and pW +> (1,3,3,4), V and W are not 
conformable. (Why?) 


Up until now, we were very careful to write expressions like 
(fT=~ppA) as ((€ppA)pl)=~ppA, so that the conformity requirements 
were trivially satisfied. For the remainder of this tutorial, we 


can use the conformity conventions just described. 
7.3 Specification Expressions 


If A is a variable and you want to store the value of an 
expression E into A, you use the specification operator <« and 
write 


A<+E 


By definition, the shape of A will equal the shape of E, i.e., 
pA «> pl. Since + is a dyadic operator, you may use it anywhere 
you would use any other dyadic operator. The value of A<E is 


the new value of 4A. 


ie UE Um aE eintivate cnet ee ee ee a a ee na et ee cee we aume wee oe a we we me ek ee ee ee aD wet wee fe 


Let E be any well-formed, array-valued expression. Then, Fis a 


select expression on E if it is a well-formed expression 


‘GEM CRED SERED GRY AEP NERD eee 
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consisting of an arbitrary number of the following Operators 
applied to E. 


(1) take 
(2) drop 
(3) reversal 
(4) rotate 
(5) subscripting 
(6) 4 
(7) A 
(8) compression 
(9) expansion 
(10) mesh 
(11) mask 


(12) transposition 


Then, if F is a select expression on the previously specified 
array A, when you write (F)+*«(E). If (eF)>pE the specification 
is eguivalent to A+F<#£. Another way of looking at this 
concept is: ee left-hand side of a specification may be any 
selection expression on A that could have been written as a 


Subscript expression on A. 


Examples: Suppose J0#G <*> 1 and A has been specified as 
A «+ (2,3)p16- Then 
ALL 2. S38). 2.412 
produces 


tf 2 4 
A «> H S 2 


dst 1: x 


1 September 1973 System Development Corporation 
| | TM-5074/100700 


This could also have been written as either 


((0,0,1)/4A) « 12 
or ( 144) « 12 


We also can write 
(QA) + (3,2)p64+16 


producing 
7 So 
A <> 8 10 12 


7.4 Conditional Statements 


APPLE, like many programming languages, contains conditional (or 
"IF") statements. A single conditional statement always 
contains some test (i.e., an expression that evaluates to a 
logical scalar). If the test is satisfied (evaluates to 1), 
then the expression associated with the conditional statement is 
evaluated. If the test fails (evaluates to 0), then control is 
iecwstouned to the statement immediately following the 


conditional statement. 


A more intricate conditional statement consists of a test with 
its associated expression and an ordered sequence of 
alternatives. If the test succeeds, its associated expression 
is evaluated and control is transferred to the statement 
immediately following the entire conditional statement. If the 
test fails, each of the alternative tests is executed until 
either some one of them is satisfied or they all fail. As soon 


as the first alternative test is satisfied, its associated 
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expression is evaluated and control is transferred to the 


statement immediately following the entire conditional 


statement. 


7.4.1 One-Line Conditional Statements 
The sinplest kind of conditional statement fits on one line. It 


is of the form 


Ik test THI] expression 


There is also a version with an alternative. This takes the 


form 


Ji @xpression 1 ELSE expression 2 


In the conditional statement, if the test is true, then 
expression 1 is evaluated and control transfers to the statement 
immediately following. If the test is false, expression 2 is 
evaluated and control is transferred 6 the statement 


immediately following. 


Examples: 
ITF A>10 THEN X « +/Y 


If 1 «>» A>icthen X is redefined as +t/YJ,. 
Otherwise xX retains its original value. 
iP A P10 PHEN XS 472 Bese fe 1 fe 


This time ¥ will be changed regardless of the value of 4710, 


se 1919: a 
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The statement is equivalent to X <« (A>10)/([/Y3+/Y) 


7.4.2 Multi-line Conditional Statements 


If a conditional statement will not fit on one line, it is 
necessary to use a multi-line conditional statement. This form 
is far more sowerrul than the one-line conditional statement, 
which cannot control the conditional execution of a set of 


statements, 


The simplest type of multi-line conditional statement involves 
Only one test and has no alternatives. It is of the form 
IF test 


statement 1 
statement 2 


statement n 
WNDLE 


Here, if the test is true, statement 1,..., statement n are all 
executed in order. If the test is false, control is passed to 


the statement immediately following the associated ENDIF, 


There is also a multi-line conditional statement with an 


alternative. 
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It is of the form 


Ik test 
s1 
s2 
sn 

BES t1 
t2 
tm 

ENDIF 


If the test is true, statements s1,...,Ssn are executed and 
control transfers to the statement immediately following the 
associated fDi‘. If the test is false, then statements 


ti,eeeygtm are executed and control is transferred to the 


statement immediately following the 44D, 


The most general type of multi-line conditional statement allows 


you to write as many conditional alternatives as you need and an 


ELSE alternative if you want one. 
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It is of the form 


wR test 1 


IS 
Peay 

IK 

thy 
(Tt 
o 
i ' 
ct 
N) 


iD 

=e 

jm 
tt 
@ 
th 
ct 
J 


But yf ? 1 
! : VW 

hoa du a g 

s eiaenatinimanatti aamectiiensanal 


EDIE 


First, if test 1is true, then statements s1,...,Sn are 
evaluated, and control is passed to the statement immediately 
following the associated #NVDIr. If test 1 is false, test 2 is 
evaluated. If test 2 is true, statements t1,...,tm are 
evaluated and control is passed to the statement immediately 
following the associated #Hi/DIF. This process continues until 
either a test is true and its associated statements are 


executed, or until all of the alternatives have been exhausted. 


If the last alternative is an £LSI, its associated statements 


= OZ, = 
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will be executed if all the preceding tests have failed. 


You should note that any of the statements associated with a 


test or an #LSE can, itself, be a conditional statement. If 


Such a conditional statement is executed, it is treated exactly 


as if it was a conditional statement occuring elsewhere in a 
program: it will either be a one-line conditional statement, 


there will be an HVDJF associated with it. 


ona Sate mnie wanton samme 


Example: 


Li (A<B)V CoD <« QR 


X«+/¥ 
Z«QA 
OR IF (A>B)aDz0 
X « AiD 
If @>kR 
2 «+ V&R 
ELSE Z + (OV)QR 
ENDIF 
ELSE 2 + (p92) 0 
ENDIF 


In this example, ) is specified in the first test. 
Regardless of the truth of that test, D will 
retain the value O+R until it is respecified 
elsewhere. The Q& I# alternative, which will be 
executed only if the first test fails, contains a 
conditional statement of its own. That statement 


determines the value of Z If both tests fail, 2 


is respecified as an array of zeros. 
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You should indent conditional statements the way we have 
done in this example. When you do, it is very easy to 
identify the statements associated with each alternative. 
This is especially true when a conditional statement 


contains other conditional statements nested within it. 


7.5 Conditional Expressions 


The right-hand side of any specification may be a conditional 
expression. A conditional expression looks exactly like a 
conditional statement, except that there is now a value 


associated with it. This is because every APPLE statement is an 


expression the value of which is discarded. 


When a conditional statement becomes a conditional expression, 
its value is the value of the last expression in the alternative 
that is executed. If it is possible that none of the 
alternatives be executed, you must be sure to include an 7,5 
alternative, otherwise the value of the conditional expression 


will be undefined. 


The shape of the value of a conditional expression is determined 
by the shape of the last expression in each of its alternatives. 
This shape is determined exactly the way it would be determined 
if these expressions were to be operated on by some dyadic 


scalar operator. That is, these expressions must be pair-wise 
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conformable. Therefore, the shape of the result of a 


conditional expression is the conformed shape of the last 


expression of each alternative of the conditional expression. 
Example: You can calculate the value of !/, 


where WiS a nonnegative scalar integer by writing 


FACTORIAL « IF W=0 THEN 1 ELSE x/1W+1-IORG 


7.6 Iteration Statements 


Iteration statements are akin to the DO-loops of FORTRAN. They 


provide a convenient means of performing the same set of 


Calculations repeatedly on some set of elements. 


where ; is a variable name, V’ is an array-valued expression, and 
S1,S52,..e-,Sn are statements. Statements si,...,Sn will be 


executed together (pV)L[ZLO#G] times, each time with / assuming 
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one of the values 


VELIORGIZORG) VIC IORG11+ZORG],..., VECLORGIT/i(pV)CLoRGI) 


exten sien wean Ree wae wen nea vane beniia wane iene qonee 


Thus, if JORG «+ 1 and you write 


DQ VE O*15 
ACI] « 2xI 
REDEA! 


{fr 


Then, this is equivalent to your having written 


AL6] « 12 
ALI23 24 
AC18] 36 
AL 24] 48 
AL30] 60 


A a + 4 


Of course, you could have writtenA[6x15] <« 2xI in this case, 
but that is because the body of the loop only contained one 
statement. Since the body of the loop can contain arbitrarily 


many statements of arbitrary complexity, you could write, e.g., 


DQ Je ips 
PF + Q(¢(2=19908)/(03;1)), (L=1ppB)/ (037) 
(FAZ) «+ AL3/(L4#1p90B)/JIOFAB 


tee awe Gee mae 


In this code, ,;1s an arbitrary array, and ,J successively takes 
on the value of each row of iopR If the dimension of B is 
unknown at coding time, this loop could not be written as a 


sequence of statements without use of some form of iteration. 
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a we ap cum Gem ae —_ ie 


sn 
REPEAL 


The test is any expression that evaluates to a logical scalar, 


and the sl,...,sn are statements. 


The test 1s evaluated and, if true, statements s1,...,sn are 
evaluated. Then, the test is reevaluated. If it is true, 
statements s1,...,Sn are evaluated again. This process 


continues until the test is false, at which time control is 


The WHILE statement is useful for controlling some process that 
must iterate while some condition is satisfied, e.g., a 
numerical approximation process. Note that iteration continues 
as long as the test is true. Consequently, you must provide a 
means for either causing the test to eventually evaluate to 
false, or make use of one of the operators described in Section 


7.6.4. 


Example: 


WHILE A/EPSILONS|X-Y 
X « ([X-16X)+2 0 | 
Y<« (Cf Y+C (OY) +V( 10Y))23 


eth abit ently amie 
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will be iterated until every component of the absolute 


value of yY-Y is less than FPSILON. 
7.6.3 UNTIL Statements 


UNTIL statements are written 


nes canine anil aime vanes 


UNTIL test 
s1 
s2 


sn 


where the test is any expression evaluating to a logical scalar, 


and s1,...,sn are statements. — 


The sequence of statements is iterated until the test is true. 
WIILE~?T . Two kinds of conditional loops are provided to permit 
programmers who think in terms of a termination condition, 
rather than in terms of a continuation condition, to directly 


translate their thought process into APPLE code. 
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7.6.4 J[teration Control Operators 


It is often necessary to prematurely terminate one or more 


iterations of a loop because some boolean condition is 
Satisfied. APPLE provides you with three anadic operators, 
[i.e., having no operands] with which you can direct the flow of 


control within such loops. 


The first of these operators is CYCLE, which 
terminates execution of the current iteration and 
transfers control to the top of the loop for the 


next iteration. 


The LHAVE operator causes control to be 
transferred to the statement immediately following 
the REPEAT associated with the loop. This causes 


00m uel wee wees ieee ee 


the loop to definitively stop iterating. 


The exif operator causes control to be transferred 
associated with the outermost loop in which the 


instruction occurs. 


either operator may be used. But, if FxiI7 is encountered ina 
nested body of loops, then all of the loops stop iterating 


immediately. 
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7.7. The CASE Statement 


The CASE statement is a simplification of the conditional 
statement. It is used when you want to evaaluate one of a set 
of Seprsasion: e6quences based on the value of some scalar 
expression, usually integer-valued. 
The CASE statement takes the form 

CASE scalar expression 


value list 1 > statement sequence 1 
value list 2 > statement sequence 2 


value list n > statement sequence n 
ELSE statement sequence 
ENLC 


madichy kien winteams drainer 


The #LSE Clause is optional. If together, the vaiue lists 
exhaust the possible values of the scalar expression, you do not 


have to provide an LSE clause. 


nem ene ene aoe 


The value list is a sequence of scalar constants, separated from 
each other with semicolons. The statement sequences may contain 
any combination of APPLE statements, including conditional or 


jteration statements. 


Example. This CASE statement computes :” and stores it in“%. 


N and Zare assumed to be one-element, nonnegative 


scalar integers. 


2 V0 = 
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CASH LOR 
O > FF © Kx/1+14 
1 > 4 © Xfi 


We could also write this code sequence using two Cf4Ssz 


statements. 


0 > Z2 «+ X/14i04 
Z<« X/iW 


Of course, this could have been more concisely 


oy 


written as Z <« x/fORCGtiN, but our intent was only 


to illustrate the use of the CAS# statement in an 


elementary context. 
7.8 The CASE Expression 


The CASE statement is a Simplification of the conditional 

statement, so any CASE statement can be rewritten as a 

conditional statement. There is a CAS! expression in the APPLE 
Qe 


language. The same conformability conventions apply to G48 


expressions as apply to conditional expressions. 


2. 
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Example: We could rewrite the second example of Section 7.7 as: 


Z+« CASE WN 
Oe ec] 
HLSE CASE LORG 
O > xK/1+i1l 
1 > “/iW 


END 
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CHAPTER EIGHT 


ARRAY OPERATIONS 


So far, we have investigated scalar operations on conformable 
arrays, and manipulative operations that produce subarray and 
permutations of the elements of arrays. Your background is now 
sufficiently strong in APPLE so that we can consider the class 


of operators that perform numerical manipulations on arrays. 


8.1 The Index of an Array Within an Array 


Suppose A is some vector or a one-component array and 4is an 
arbitrary array. Then 41%, the index in4A of 4 is an array such 
that 


SHAD. > 06 


and for each row Lof 1ep2, (Ai2)L;/L2) is the least index / such 
that (,A)C7] «> BOL;/L)]. If bL;/L2] is not an element of 4, then 


(AiB)CL3/L2] «> 14+f/i9,4a . 
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Examples If JOU «+? 1, A «> (3, 7,2,0,1,4%,1) and 


41 3 2 
“J 4 5 
iio«-> 
Oo 41 4 
2 
Then 
5 1 3 
2 6 8 
AiB +> 
uoe5 6 
1 8 5 


Note that 4[5] «+ 4A[7] <* 1, but the index 
returned to 418 is always 5 since 5 is the least 
index J for which ALJ] +> 1, Note also that an 
index of 8 was returned for the two elements 5 and 


"2, which are not contained in 4. 


You can tell if every element of / is an element of A since, if 


that is the case, we must have 


1 «> v/(,A1B)SI/1p9A 
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8.2 Array Membership 


Let A and B be two arbitrary arrays. Then, 4€8 is a logical 


array of the same rank as A and contains a 1 corresponding to 


each element of 4 that is present somewhere in 2. 


Example: Suppose 4 and » are the arrays from the example in 


Section 8.1. 


Then 
Ac& «> (1,1,1,1,1,1.,1) 
and 
+. 1 x 
aw a) 
BeA +> 
dk, A 
de “> St 


8.3 Sorting 


If Ais a vector, then AA is a vector, p4A <> pA, Such that 
(AA)[ r]tS the index of the J-th smallest element of A. That is, 
ALAA ]LS a vector whose first element is the least element of A 
and each of whose remaining elements is no less than its 


predecessor element. 


Just as 4A can be used to sort the elements of A into ascending 


order, YA +> 4A can be used to sort A into descending order. 
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Examples: If JORG «7 tand A «>(6,3,5,3,9,3,1) 


Then 
AA «> (7,2,4,6,3,1,5) 


¥A +> (5,1,3,6,4,2,7) 
8.4 Quter Products 


Suppose Aand sare any two arrays and #®is any scalar dyadic 


operator. The »x outer product of 4Aand B, written Ac.*f, is an 


RED: GED NGS KER EH BS! EP SUE EE PE CELE 


array containing the x product between every element of A and 


every element of 2B. 


Formally, pAc.#B «+ (pA),p2 and for each row L of 1pA°.#d, 


(Ao.#B)C3/L] <> AL3/(ppA)tLI#BLs/(-p 0B) tL) 
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Examples: Suppose 4 «+ (1,2,3,4,5,6,7) 


Then 
1 2 S mt 5 6 9 
2 4 6 8 10 12 14 
3 6 9 12 15 18 21 


ie en ee i SS AG: OO Oe Os 


41 0 0 0 0 6) O 
O 41 0 0 @) 0 0 
0 0 ah Q @) 0 @) 
Ao,=A ~«> 0 Q @) 1 0 ) 0 
Q O 0 0 1 0 O 
0 6) ¢) @) 6) a: O 
Q 0 O Q 0) @) 1 
i Fi i 1 41 1 41 
0 si 1 gf i) 1 
@) 0 1 8 r al 1 
Ae.sA «+> 0 Q 0) 1 1 yi 1. 
O Q Q 0 1 4. dk. 
9) 0 @) 6) 6) 1 1 
0) 0 0 0 6) 0 l 
@) 6 6) 
OQ? As2) GZ 
O AS 23 
C(253)606512,45585. 5,2))°.0(0,67513) <s 
0 8 8 
“5 °§ 5 
0 2 2 


8.5 Inner Products 


The APPLE inner product is a generalization of the linear 


algebra inner product of two matrices. Suppose that 4 and B are 
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two matrices such that 1+pA «> 1toB, The following code 
produces the inner product Cc of A and 2B. 
C+#((1+9A), 1tpB)p0 
DO L € ipc | 
COh3/L)<«+/AC14+23; ]*80 3145) 
REPEAT 
(In normal parlance, ClI;:/] is the sum of the componentwise 


product of the I-th row of A with the J-th column of 8B.) 
The APPLE notation for this inner product is A+.~xB. 


In general, if 4 and 8 are matrices satisfying 1+tpA <> 1+p8 anda 
* and *® are two scalar dyadic operators from Table III (see 


Section 5.3), then the *-* inner product of A and 8 is written 


Ax.MB , where pA®.aB «> (1tpA), 1tpB and for any row L of 1pAx. «2 
A®.MBDL3/L] <% ®/ACL1tZ3;]«BL 3145) 
Examples: 
1 2 - Bg 9 422 -2o 
3.4 +.*% 4 5 6 +> 19 26 33 
5 6 29 40 Si 
ds, 2 e238 oS Ff 


OO 
Bs 


+.F§ 45 6 +* 7 86 @Q 


In the remainder of this section we generalize the inner product 


to conformable arrays of arbitrary rank. 
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If ppA «> ppb «> 0 then Ax.KB +> &/AMKB. 


Otherwise, Ax®.“B is defined only ifA and B&B satisfy one of the 


following conditions: 


(i) pp0A 
(2) — ppB 
(3) “AtpaA 
(4) 1+oB 
(S) “1to“ 


8.5.2 Definition of 


If the conformability 


conceptually replaced 


if .6oA 

A 

Li 

(2°) Tf pp 
A 

B 


(3") If “1tpoA 


ICs bs 


(4") If 1tp8 


JOS [ts 


(5') If 1+pA 


the General Inner Product 


conditions are satisfied, then A and 


by 


<> 
“> 


<-> 


oo. 


<--> 


<> 


4-> 


<> 


arrays 4 and B given by: 


~t4+pA)pB 


1te8 then 


A 


, 
1 then 


A 
(( itpA),14pB)pB 


then 
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A <> (d1pp4 )Q((Qt0B),> 149A) p(bippd )&A 


B «<>? &B 


Then Ax.«E +> AxX.KB +> C where 


pC «> ( 1494),140B 


and for each row Lof 10¢ 


Cl3/L] «+ */(GAA)«HAB 


where 
CLsl/i2) +> CC t4ppA)o1)50 
Clal{i2) «UC Lappe) tL), 0 
Wisb/12] <+ 0,¢€ 1+ppB)p1 
HE;sf/12] «> O,(Ci-poB) tL 
Example: 
“3. °°5 
e/a: 
1-2 1-2 a ee 
Came 3 4 3 °° 5 
oS <--> 
Hy 6 5 6 i) _3 
7 8 7 4 


mak 
VAe 


8.6 Change of Base 
The number 2307 = 2x103+3x102+0x10!+7x109. This is a polynomial 


representation of the number 2307 in the radix 10. We also have 


2307=3x5%+3x534+2x5241x51+2x59 as a polynomial representation of 
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2307 in the radix 5. 


If you have a vector representation YVYof a number WV in the base 
B, yOu can compute the value of J by evaluating the expression 
(x\(( 1+pV)pB),1)+.xV. For example, if B <«>5 and 
V «++ (3,3,2,1,2), the expression evaluates to 
(x%\ 05,555,551) )46%(39435251,2) <> (6255125, 2555.1)4.*03, 3424152 
“> 18754¢3759+504+7 <> 2307 
The APPLE notation for this conversion is B.iV (read "8B decode 


Very , 


Conversely, if you want the vector representation V of W to the 
base 8, then you can compute JV by the Euclidian Algorithm as follows: 


Vee 
UNTIL N=0 
V+(BI|N),V 


The APPLE notation for this inverse conversion is BT (read "B 


encode WV"), 


Change of base can be generalized to cover mixed bases. If B is 
a vector, then 281V is defined when B and V are conformable for 
inner product: 


BaV +> (x\(14B),1)+.xV 


Indeed, Bi1V is similarly defined for any arrays B and y that are 
conformable under inner product. The formal definition is 
extremely complicated, so let us look, instead, at a few 


illustrative examples. 
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Examples: Suppose you want to know the number of seconds 
in 1 year, 20 days, 5 hours, 3 minutes and 17 
seconds. We will proceed with the assumption that 
1 year=365 days, 1 day=24 hours, 1 hour=60 
minutes, and 1 minute=60 seconds. Then we would 
have to compute | 
(365x24x60x60),(24x60x60),(60x60),60,1)+.x(1,20,5,3,17). 
But this is equivalent to 


(x\(365,24,60,60,1))+.*x(1,20,5,3,17) 
+> (1,365,24,60,60)T(1,20,5,3,17) «> 33282197. 


Note that the first element of the vector (1,365,24,60,60) could 
have been any arbitrary number, since it is 
discarded in the evaluation of 1. 


Similarly, (1,3,12)1(3,1.6) +> 126, the number of 
inches in 3 yards, 1 foot, 6 inches. 


Neither B nor V need be integral. 


24 60 60 763 430.47 110 22563 14545.47 3660 
1 3 12 2 208 118.97 30 <> 327 179.07 48 
1 20 12 43 28.07 6 1755 995.07 252 


The inverse operation BtW is also defined for arrays B and jy. 
The result is an array VY such that Biy <> W. Hence, 
oBTN <> (9B8),pN and the “base Bg vectors" run along the first 


coordinate of the result. 


Examples: 
(505)T24U42 <> (3,4,2,3,2) 
(59 2) T 13 €> (1,1,1,0,1) 
(192) T 13: <> 00,0,0,.1,1,0,1) 


(1780,3,12) T 126 «> (3,1,6) 


This last result is the number of yards, feet and 


inches in 126 inches. 
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ee OD 


(10,10,10,10) T (456,7890,2345,123) «<> 


NO 


0) 
m 
5 
6 


Ow won 
im £— W ho 


- 
Me 
‘ae 


8.7 Matrix Inverse Operator 


Suppose B iS a nOn-Singular matrix and 2/pB +> 1 So that 4 has 
at least aS many rows as columns. Then, $2 «> £ such that £4+.xs 
is the identity matrix. i is the left inverse matrix of & and 


pli <s: OOB. 
8.8 Matrix Division 


When you have a vector or matrix 4 and a matrix % such that 
2/oBb +> 1, then AHS is defined when itpA +> 1itp8. By definition 


Example: 
109: “2 uy 8 


5 » 
97 56 §} 3 9 2 «++ 9 3 
114 87 a 40. 2 5 4 


= 123. 
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CHAPTER NINE 


SUBROUTINES, FUNCTIONS AND OPERATORS 


This chapter contains a discugsion of how to write subroutines 
and functions. The difference between the role of functions in 
APPLE and those in other programming languages is that there is 
little conceptual difference between functions and operators. 


consequently, APPLE can be treated as an extensible language. 


9.1 The Distinction Between Functions and Subroutines 


A subroutine is a code sequence that is parametrically self- 
contained. It can be invoked from SHY park of a program. After 
the subroutine has completed execution, control is returned to 
the first executable statement following the point from which it 


was called. 


A subroutine may use or manipulate the contents of variables 
used by the main program, Or it may use or manipulate variables 
that are accessible only by the subroutine itself. If a 
subroutine is parametrized, then the values for the parameters 


are specified at the points of the program at which the 
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subroutine is called. 


A function, on the other hand, always has a value associated 
with it. Like a subroutine, it may use or manipulate variables 
used by the main program, or it may use or manipulate variables 
accessible only by the function itself. A variable can be 
assigned the value of a function. But, in order for a 
subroutine to modify variables external to itself, it must 


explicitly assign values to them. 


Subroutines and functions may be anadic, monadic, dyadic or n- 
adic in that they may take no arguments, one argument, two 


arguments, Or n-arguments, respectively. 


If yorrcekK is an anadic function or subroutine, you invoke it by 
writing yorrck « If pozonrysgs 1s a monadic function or 
subroutine, you invoke it by writing POLONIUS A, where A is some 
expression representing the value of its parameter. If 
ROSENCRANTZ is a dyadic function or subroutine, you invoke it by 
writing A ROSENCPANT? B, where aA and B are expressions 
representing the values of its two parameters. If GUILDENSTERI 
is an n-adic function or subroutine, you invoke it by writing 
GUILDENSTERI ¢ AsB3C3D3...63;Z2>, where A,B,...,Z2 are expressions 


representing the values of its n parameters. 


Note that the monadic and dyadic functions and subroutines are 
written in precisely the same way monadic and dyadic operators 


are written. Since functions have values, they can be invoked 


~ 125 - 


1 September 1973 System Development Corporation 
TM-5 0747100700 


from within any expression. Consequently, there is no 
difference between the use of the functions you define and the 
APPLE operators. In the remainder of this chapter, we will 


discuss how you may define and use your own operators in APPLE. 


9.2 The Form of a Function or Subroutine Definition 


The definition of a function or subroutine consists of a 


heading, a body, and a footing. 


The heading names the function or subroutine, specifies whether 
it is anadic, monadic, dyadic or n-adic, whether it is a 
function or a subroutine, and names all of the variables local 
to its body. The body consists of the APPLE statements that 
perform the computations and manipulations of the function or 
Subroutine. The footing ideminaces the definition of the 


function or subroutine. 
9.2.1 The Heading 


The heading begins with the symbol Y. If you are writing a 
function, you subsequently specify the name by which you will 
refer to the result in the body of the function, then the 


specification symbol <. 


If the function or subroutine is to be dyadic, you next list the 


local name of its left argument. 
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The following entry is the name of the function or subroutine. 


Then, if the function or subroutine is monadic or dyadic, you 
list the local name of its right argument. If it is n-adic, 
then you enclose the local names of the parameters in the 


brackets Cand >, separating the parameters with semicolons. 


Finally, if there are to be any variables local to the body of 
the subroutine or function, you list their names, preceding each 


one with a semicolon. 


Examples: The heading for the anadic subroutine YOAICK is 


V LOR LCA 


This subroutine has no local variables. The 
heading for the anadic function HORATIO, which has 


two local variables Mand J, is 
Vo in & AORATLO sii $i) 


Here, the value HORATIO returns is to be explicitly stored 


in a local variable called 2. 


The heading for the monadic function POLONIUS is 


V2 ¢ POLONIUS Y 


Here, Y is the local name of the right argument of PoLONIUS 


and the value of POLONIUS will be stored in 2Z. 


If you wanted to define a monadic subroutine OPHINLIA, you 
would write 


VY OPH TIA 
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OPHELIA was not defined as having any local variables. 


A dyadic function ROSENCRANTZ, with local variables 


6 oe are eh would have the following heading 
Vg e& XK POSTUCAANT 2. J eQ sneer 


Here X and y are the left and right arguments of 
ROSENCRANTZ , the value of which will be stored in the local 


variable 7. 


If GUILDENSTERN is to be a S-adic function of 4,7?,C,D2,E and 
has a local variable 7, you would write the following 
heading: 


V4 « GUILDENSTERN CA3B3C3D3;#23N 


9.2.2 The Body 

The body is an APPLE program. If the subroutine is a function, 
the result must be stored into the name you specified someplace 
in the body of the subroutine and in the heading. It, and any 
other variable names from the heading may be used as local 


variables within the body. 


The body may also use the names of variables existing outside 
the body, i.e., they are neither the result, operands, 
paraneters, nor declared local names. Such variables may be 


modified from within the body. 
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9.2.3 fhe Footing 


The subroutine of the function definition is terminated with the 


symbol y on the line following the last line of the body. 


9.3 Call by Value 


The parameters of a function or subroutine are evaluated from 
right to left before it is entered. The resulting values are 
copied into the local storage area of the routine. Then, the 


body is evaluated. 


Therefore, it is impossible to modify a variable existing 
outside the routine by using the variable as a parameter and 
then trying to modify it inside the body of the routine. The 
parameters must be considered as local variables that have been 


initialized when the routine was called. 


9.4 The Scope of Names 


Inside a body, you may manipulate global variables defined 
inside other routines. If several variables have the same name, 
Only one of these is accessible. In order to determine which 
variables are accessible, you follow the chain of calls 
backwards. The first time a specific name is encountered, you 


have found the one that is accessible. 
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Suppose that A,B and C are names of global variables. 


There are defined routines Ff1,F2 and F3 with the 


following headings: V4<¢« Fi Aj D Suppose 
7 Vv. A F2 ¢ 
V A+ F3 D 


that Fi is called first. Then, inside the body of 
Fi, variables with the names 2Z.4,2,C and D are 
accessible. 8 and C are global; Z and D are 


local, as is A since it is the parameter of Fi. 


Assume that F2 is called from within the body of Fi. Then, 
the body of F2 may operate upon its own variables A and C, 
the local variables Z and Din F1, and the global variable 


Be 


If the body of F2 calls F3, the body of F3 may operate upon 
its own local variables A and C from F23; the local variable 


Z from Fi; and the global variable B. 


If the calling Sequence was different, quite new effects 
might appear. For instance, if F3 is called first, it does 
not have access to any local variable Ze If the body tries 
to use Zin an expression, an error will result. The 
variable A is no longer local to F1, but is the global 


variable by the same name. 


9.5 Recursion 


Functions may call themselves recursively. In such cases, each 
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incarnation of the function may be considered as a separate 


function with its own environment. In particular, the 
Function's local variables are local to it and are global only 


to those routines it calls non-recursively. 
9.6 The RETURN Operator 


If yOu want to exit from a function prior to its completion, you 
may use the pfrryRy Operator. This operator may be used either 


anadically or monadically. 


ome GE Cm ame gua om 


specification expression, which sets the return parameter. For 
example, you could write: 


LF A # 0 THEN RETURN Z< ?A 


9.7 Comments 


Any line in a routine's body can be made into a comment if you 
start it with the symbol a. The comment symbol is called the 
"lantern" because it often sheds light on a complicated 


program. 
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9.8 Implied Loops 


Suppose you define a monadic functionOSRIC . When you write 


OSRIC A, the function is applied to the array Ain its entirity. 


If you would like to apply OSRIc along coordinate 7; of A, you 


write Z+<+ OSRIC [I] A. 


Then, OSRIC will operate on the hyperplanes along coordinate J 


of A, producing a resultan+ array such that 


pOSRIC [I] A <> ((I-IORG)t0A),(pZ),(I+1-LORG)+pA 


where for rows J of i(pA)LTI] 
ZCLIORGI]J] +> OSRIC FAA 


and 
F «> Q((T=1ppA)/(031)),(LT=pp4)/(0327) 
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CHAPTER TEN 


ON THE ORDER OF EVALUATION 


The order of evaluation in APPLE expressions is controlled by 
the rules of right association. However, you can easily 
construct examples in which it is undesirable to strictly adhere 
to this philosophy. For example, if we were to write 


Ko C1e0,1 4050) TERATELYCOMPLICAT  DEAPRE OOLON 


we see that X will receive only 40 percent of the value that was 
computed on the right side of the specification symbol. If 
there is no specification of any part of the remaining 60 
percent of TERRIBLYCOMPLICATEDEXPRESSION , then we have wasted 


considerable computing power in its evaluation. 


Consequently, there is a slight, but important, modification to 


the principle of right association in APPLE: 


No portion of a select expression is evaluated unless 


it is subsequently stored into a variable. 


This principle has an interesting consequence. Consider the 
code: 

B £> C150 e2.0) | 

Y«<« (B#0)\CB40)/4+B 
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Now strict adherance to right association yields disastrous 
results since the machine must Givide 4 by 0 twice. But 
suppose, instead of producing an interrupt, the result of 
division by 0 is the undefined scalar [{. Then, 

ueBe> (4,(11,2,0)) . Continuing the evaluation, we see that 
(B40)/(4,0,2,0) «> (4,2). At this point, the unwanted [) has 
disappeared. Next, (B#0)\(4,2) <> (4,0,2,0), Thus, when all 


has been computed, we have Y <> (4,0,2,0). 


This was undoubtedly what the programmer had intended. APPLE's 
modification to the principle of right association tends to 


correspond to what programmers find natural. 


The APPLE compiler maps all selection expressions into a 
standard form (given in Section 1.11 of the APPLE 
specification). The compiler is then able to distinguish the 


necessary and unnecessary computations, suppressing the latter. 
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APPENDIX I 


AN INDEX OF SYMBOLS 


This appendix is an index of APPLE operators. 
It is divided into two parts: arithmetic 


operators and array operators. 


MOWADIC 


Symbol 


Name 


identity 
negation 


siqnum 


reciprocal 


exponen- 
tial 


natural 


logarithm 


£Lloor 
ceiling 


absolute 
value. 


roll 
negation 


factorial 


pi times 


ARITHMETIC OPERATORS 


Page 


18 
19 


19 


20 


20 


23 


IA 


IV 


DYADIC 


Symbol 


A=? 


Name 


addition 
subtraction 


multipli- 
cation 


division 


exponen- 
tiation 


loqarithm 


maximum 
minimum 


residue 


deal 
combination 


circular 
function 


AND 

OR 

not equal 
equal 
NAND 

NOR 

less than 


less than 
or equal 


greater than 
greater than 
or equal 


25 


25 


26 


26 


28 


29 


27 
27 


26 


35 


34 


29 
30 
31 
31 
32 
32 
32 


33 


33 
3 4 


33 


MONADIC 


id 


Ce erailaonthemmalornnnkamaed 
ernie comme cemre Manel tee 
Vanes came wane GE ume 


Name 


iota 
odometer 
ravel 
shape 


partial 
subscript 


reduction 
accumulation 
interval 
transpose 
reversal 


subscript 
generator 


whole array 


matrix 
inverse 


123 


98 


101 


105 


108 


107 


109 


109 


109 
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ARRAY OPERATORS 


DYADIC 


S ol 


in 


Name 


index of 
element of 
row of 
laminate 
catenate 
reshape 
compression 
expansion 
subarray 
cross section 
mask 

mesh 

prefix 

suffix 
transposition 
take 

drop 

rotate 
specification 
grade down 
grade up 
outer product 
inner product 
decode 

encode 


matrix divide 


Page 
113 


115 
105 
85 
64 
57 
68 
65 
65 
67 
78 
79 
72 
72 
87 
Ge, 
76 
+2 
96 
115 
115 
116 
117 
120 
121 


123 


APPENDIX II A Fast Fourier Transform 


In this appendix a simple program that performs a Fast Fourier Tranform is 
delineated and described. It is assumed that you are familiar with the way 
FFT's work, (If not, we recommend that you read "A Guided Tour to the Fast 
Fourier Transform," by G. D. Bergland, IEEE Spectrum, July, 1969, pp. 41-51. 
This paper includes a comprehensive bibliography. ) 


The algorithm we follow is an adaptation of a method developed by R.D. Schmidt 
and communicated to the author by W. Juran, Proprietary Computer Systems, Inc., 
Van Nuys, Cal. 91406. 


We assume that the input data is an a rank-2 array B, where pB < (L,2), with 
ZL an integral power of 2 not exceeding some fixed number, say, 256. Each row 
of’ B consists respectively of the real and imaginary parts of a data value. 


First, we produce a simple preprocessing function. The algorithms for trans- 
formations for time to frequency and frequency to time are essentially the 
same. The only difference is that in the time to frequency transformation, the 
result must be divided by the number of data sample points. 


With this in mind, we can write our simple driver function, FAST. We will 
make FAST dyadic: the second parameter is the data, while the first parameter 


A determines that the transform is time to frequency if A «+» 1 and frequency 
to time otherwise. 


The code for FAST is shown below, where we presume that JORG <> 1: 


V Z<« A FAST B3C 


C1] IF ~ ((C + $9B)02])«0,18 LZHEN RETURWT Z + (pB)pERR 
[2] Z + QCoFFT((1+2@CL2])p2)pB 

tod D2) eG 

[4] aFFT RETURNS CONJUGATE OF RESULT 

[5] IF 1=A THEN Z + Z:C([2] 

v 


Here, we insist in line [1] that the number of rows in B be a power of 2 less 
than or equal to 256. If not, FAST returns an array Z containing a predefined 
error value contained in ERR. | 


In line [2], the monadic function FFT is invoked, returning an array whose 
transposed value is stored in Z. We will look more closely at line [2], 
momentarily. 


Line [3] produces the complex conjugate of the result of FFT. 


And line [4] divides by the number of data points if the transform was time 
to frequency. 


When FFT is invoked on line [2], its argument is a rank 1+2eC[2] Pend 
of B such that ,B[[1]1] consists of the pure real data components and 

»BLL1]2] consists of the pure imaginary data components. For example, if we 
had started with B <> Q(2,16)(116),-116 then the argument X transmitted to 

FFT would have been 


X «> (2,2,2,2,2)p&B 
1cC oy - 


ALISA 5s 4. eS 


The imaginary components X[2;2;;;1] and X[2;2;;;] are just the respective 
negatives of these. The motivation for this restructuring will soon become 
apparent. 


The FFT must take care of the required binary sortings on this array. The 
first such sorting is the one in which the real and imaginary components 
respectively assume the position determined by reversing the binary encoding 
of their index in ,X[[1]1] and ,xXLL1[2]]. This can easily be achieved by 
the transposition, since each coordinate of X is indexed by either 1 or 2. 
The general desired array Y is giwen by Y <> (1,61+11 2@x/1+4pX )&X. 


In terms of our example array X, this is 


Y<> (1,5,4,3,2, )QX 


Vid ee a, oS 


cD 
— 
oy 
wi 


Again, the imaginary component is symetric to the real component. 


Next, we need the appropriate array of coSines and Sines for the real and 
imaginary components. These are given by the array 7 where 


ey (2,1)°, 000 O1M-1)QC(CM-1)p2)p0,(1 147 2)* 22H 
and 
Mo <> L2ex/1+pX 


TY o<> x/1+¥pAa 


That is TLL1J1] consists of the cosines of the appropriately transposed 
array of multiples of 2n divided by the number of data points; i.e., 
0,(02#N),(O47N),...,(02x 14+W¥+2)+N. TLL[1]2] consists of the sines of the 
appropriately transposed array of the same multiples of 02:N. 


Before we proceed any further, we list the function FFT. 


VY « FFT X3d3Q3N3N37 


[i] M<«| 2@V<x/1+0X 
[2] Y<(1,01+1M) 0X 
ES: T+(2,1)°.,00(O1M-1)Q((M-1)92)p0,(1 14+NF2)x22M 
[4] J+2+ppLr 

WHILE 1 
[6] Q+«+-/[Jy<J-1] Y 

Ye(+/(LJ] Y):€7] (-/£1] @x@):01] +/£1] 7xo[1] Q 

[8] LF JsitpY THEN RETURN 
[9] Te(1, 164+1M¢-1)07,02] 7+(1,0)/C2] 7 
ETO REPEAT 


Lines [1] - [3] contain the code for defining the initial sorted arrays 
Y and fT. 


The loop spanning lines [5] - [10] iterates until the condition on line [8] 
is satisfied. The variable J is used to control the processing of the data, 
starting with the last coordinate of Y and concluding when the second coord- 
inate has been processed. Y and fT are restructured in lines [7] and [9] 

to conrespond to the sortings required by the algorithm. 


The easiest way for you to understand the way the algorithm works is for you 
to try following its execution on a small array, say one with eight elements. 
You can then use a form of mathematical induction to verify that it works on 
the higher~-dimensional cases. . 
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Individual differences and group interaction. Since each evaluator will 


have his own set of values (S,'s)> there may be collected diverse values of 
S, for the group of evaluators. This divergence can come from differences 

in RY values and/or w (weights) values. The Delphi technique or simpler 

group interaction may be used at this time. The Delphi technique may have 
been used earlier for both RY and w values, but the evaluators are not likely 
to come to a complete agreement on one set of values (but eventually there 
will be fewer sets than the number of evaluators). Any group interaction may 
yield some influence toward agreement. Some studies on group planning seem to. 
indicate that people may widely disagree on objectives and criteria but may 


readily agree to favor certain alternatives. 


Let us see how ierevent sets of 5,8 for different evaluators can be 
compared. Let 8 represent a composite attribute of all a, 'S that have been 
considered. In our park example, 8 will be "acceptable park design with all 
the features properly balanced." S, values can now be used to produce 
Re (A,) for each evaluator automatically (associating 4 values with numbers 


in the interval [0, 1]). A display of R,'s for all alternatives A> Ay» ee 


B 
A, can be made for each evaluator (see Figure 4). The total display of all 
such values of all the evaluators may be shown in a scrambled order to maintain 
anonymity, if desired. In addition, a statistical group response, such as 
quartiles (Q1> M (median), Q3) of each R,(A,) value, can be calculated. 


Seeing where his own evaluation stands within the group response may aid him 


in understanding the overall evaluation. 
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Figure 4. Overall Evaluation of Alternatives 
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As an evaluator studies ‘the relative merits of alternatives and gains a deeper 


understanding of trade-off implications, his total conception of the situation 


grows and his fuzzy-set mapping may become more defined. If Re is his 
(t) 
preference function operating on the fuzzy-set at time t with respect to the 
global criterion imposed by 8, then as t increases, R tends to converge 
(t) 


toward a more precise function; i.e., the evaluator becomes better able to 


sort out alternatives. 


Fuzziness (or impreciseness) of many decision-making situations is qeaeray 
caused by a mixture of ignorance, randomness, and intrinsic fuzziness. This 
kind of exercise can help the evaluator to separate out types of fuzziness 
involved and learn to identify where more information is needed (case of 
ignorance), where probablistic treatment is needed (randomness), and where 
increased awareness of other value systems besides his own is needed 


(intrinsic fuzziness). 


In complex decision situations where many competing factors must be properly 
accounted for simultaneously, the interactive system can be made to keep erable 
of the evaluator's tendencies. For example, suppose the evaluator is 
excessively cost oriented and his assignments of grades of membership for the 
cost attributes fall consistently outside the interquartile range (Q, Q) of 
the group response. The system can remind him of other important factors and 


trade-off considerations. 
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Interacting with other evaluators through the system and seeing where his own 
evaluation stands within the group response may influence him Sop eeke additional 
factors into account or to adjust his preference functions. If he feels 
strongly about his preferences, he can try to persuade others by stating the 
reason why the value should be lower (or higher) than the values (in the 


interquartile range) expressed by the 75 percent majority. 


The process of interaction and reevaluation can bewecuneed until, it is to be 
hoped, convergence is attained. Polarization may occur but completely flat 
distribution is not likely for most value-laden questions. Individuals seem 
more responsive to value-oriented questions than to factual questions; that 
is, changes of opinion seem more readily attainable on value-oriented issues 


than on the factual ones. 
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SUMMARY STEPS OF MACHINE-AIDED EVALUATION 

There are seven evaluation steps to be considered, some of them will be 
elaborated on in the following pages. For simplicity, the procedure ee 
a single user (evaluator) and group interaction is not emphasized. Steps are 


presented in the "usual" order but can be reordered at the user's direction. 


1. List alternatives by name or number assigned to each. 

2. List criteria for evaluation in terms of attributes. 

i ee Rank attributes and assign weights. 

4, List values in their "natural" description (numerical or non numeric) 
for each alternative's attributes. 

oe Determine grades of membership of all values of attributes. 

6. Calculate the summary value of each alternative. 


iE Repeat any or all the steps above. 


The list of alternatives and attributes may be prepared in advance 

covering steps 1, 2, and 4 and can be thought of as an attribute-alternative 
table (Figure 5). Unlike mathematical tables, this table can contain both 
numerical and nonnumerical descriptions, even lengthy discussions supple- 
mented by pictures that can be referenced. Therefore, the physical form of 


the information may not look like the table in Figure 5. 
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Alternatives ‘ 


Attributes 


Figure 5. An Attribute-Alternative Table 


Following is an expanded description of the seven evaluation steps. 


Ly List alternatives by name or number assigned to each. 


2% List criteria for evaluation in terms of attributes. 

Attributes can be given on a noncommittal trial basis with full recognition 
that they are likely to be iradequave or incomplete; or they can be carefully 
selected by a group of people (e.g., noitey makers, planners, experts, 
representatives of the public). Attributes may be separated into two groups, 
"desirable attributes" and "undesirable attributes", or they can be all mixed 


together. Subsequent instructions will reflect the choice. 


3. 


* 
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Rank attributes and assign weights. 


Attributes are rank ordered in terms of their relative importance in 
contributing to the objectives. If there are many attributes and ranking 
is difficult, machine assistance can be provided by showing only two 
attributes at a time to the evaluator. Judging the relative importance 


of two attributes is much easier than ranking the whole list.* 


If the evaluator's judgment of importance is transitive and total in 
ordering, the attributes are listed in the order of importance, possibly 
placing two or more in the same rank in case of a tie. The logic 
(computer programs) inside the machine can easily check any inconsistency 
and ask the evaluator to compare again those attributes whose rankings 
are in conflict. If the inconsistency is not removed, it is likely that 
at least one attribute should be redefined in terms of two or more other 
attributes, or some attributes should be grouped together as one. An 


experienced evaluator usually can sense which ones are in need of adjustment. 


Showing two attributes at a time makes comparison easier but the 


resulting ordered list should be reexamined as a whole to guard against 
any possible context shifts, which may result when only two items are 
compared at a time. 
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When ranking is complete, the attributes are displayed in rank order. 

The evaluator is now asked to assign weight w, to each attribute, starting 
at the top-ranked one(s) with the weight of 100. Using this as the point 
of reference, the other attributes are also assigned weights. These 
weights should reflect the relative "strengths of effects” of attributes 


contributing to the objectives. 


4. List attribute values in their "natural" description (numeric or non - 


numeric) for each alternative's attributes. 

If the attribute values are given in verbose descriptions, they can be 
condensed to a few key words to be displayed along side the names of 
attributes. The original information sheets should also be available 


to the evaluator. 


5. Determine "grades of membership" of all che attribute values. 
Using the fuzzy-set concept, each attribute value is judged in terms of 
"grade of membership" -- i.e., a number in the interval [0,1]. Since 
comparability is important in value guaeenene: other values of the same 
attribute from different alternatives can be shown one or two at a time 
for comparison. In Figure 5, this process corresponds to moving hori- 
zontally across the alternative on the — attribute line. When all 
the attribute values are judged in this way, a new table of values is | 
created within the uniform scale. When it is displayed, it will look 


exactly like Figure 5 containing a single number (between O and 1) in 
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6. 


7. 


each cell of the table. Seeing the total array of numbers may prompt him 


to change his earlier choice of values. 


Calculate the summary value of each alternative. 


The summation S, = 2 we RY (A,) is calculated for each alternative 
i=] i : 


ae as one of the basic machine aids, but some other forms of getting the 
summary value may be tried out. The evaluator can specify his own ideas 
easily with the man-machine communication language, User Adaptive 


Language (UAL) (see Hormann, et al [1970]). 


If attributes have been separated into "desirable attributes" and 
"undesirable attributes," 55's are calculated using only those a's 

in the desirable category, and 5 = w, (1 - re ae is calculated using 
only those a's in the undesirable category. Weights, attached to attri- 
butes, remain the same since they should reflect the "strength of effects" 
regardless of desirability or undesirability. The difference, ay ~ S; 

for the alternative A, may be called the "net-benefit value." Comparing 


these values presumably will indicate a tentative conclusion concerning 


which alternatives are best. 


Repeat any or all the steps above. 


The evaluator is encouraged to go back and examine his previous judgments. 
it is usually advisable, the first time around, to use first impressions 


in making attribute rankings and in making a judgment of grade of 
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membership without too much deliberation. Stepping through the whole 
sequence rather quickly the first time, rather than dwelling on a single 
factor in detail, will give him a better understanding of how certain 


factors are accounted for in the total evaluation. 


Iterating the evaluation process tends to bring many assumptions into 
the open, and the evaluator may become more aware of how the conclusions 
are related to the assumptions. For example, assumptions on the objec- 
tives will influence the interpretation of objectives and criteria in 
terms of attributes and will also influence attribute ranking and weight 
assignments. Assumptions on political and technical constraints on the 
proposed designs will certainly influence many decisions. Probing into 
them with "what if'' questions may separate out "real" constraints from 
imagined ones or those that can be overcome by negotiation or by 


creative problem solving. 


The evaluator may, in the light of new insights and understanding, wish 
to redefine objectives and specify relevant attributes more carefully. 
Interacting with the other evaluators, or even with the policy elec 
may bring further clarification. Possible use of the Delphi technique 


has been discussed, so it will not be dealt with here. 
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Complex trade-off implications, which are typically nebulous, can be 
made clearer if bar graphs such as those shown in Figure 3 are used. 
They can be rearranged to show the Ro. values of different alternatives 
horizontally for each OL « The re may be encouraged to ask "what 


if" questions on possible trade-offs that are not evident in the design 


alternatives; the answers may suggest a new or modified design. 


POTENTIAL APPLICATIONS 


There are a number of areas of potential application for this method. 


fi Complex equipment with many performance criteria. 


Evaluation of different designs of complex equipment such as aircraft 
and underwater exploration vehicles can use machine-aided evaluation. 
In these, many attributes must be included in evaluation and they cover 
both factual information and subjective value information. This class 
of problems is less fuzzy not only because factual information tends to 
dominate but because the physical boundary in measuring operational 


behavior is relatively clear. 


2 Selection of suitable locations for large complexes. 


The problem of selecting a suitable location for a large complex, such as 
a new housing development, often requires careful consideration of many 
attributes that are qualitative in nature. Among many possible locations, 


one or a few candidates are usually selected in order to proceed with 
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designing, legal and financial negotiations, etc. Other examples of 
large undertakings whose location Saeiatons tend to affect various 
segments of our society are: manufacturing plants, airports, hospitals, 
health-care eenvers, sanitoriums, rehabilitation centers, educational - 


institutions, trade centers, highways, and transportation networks. 


3. Complex combination of things that interact. 


Hakiag a appropriate EDP sabe selection from all possible combinations 
of available hardware/software products to meet the user needs is a complex 
problem. Guessing at a suitable hardware/software mix is hard enough, 

but evaluation of a wide number of configurations when the components 
interact usually requires an advanced modeling technique (Sutherland 
[1971]). Information on the performance hae octeriactcs-iot hardware and 
software components are separately available, but very little information 
can be had on the total performance characteristics for specific 
configurations--unless the pieces are all ee ee same manufacture. After 
modeling produces the system's performance characteristics, our technique 


can be used in total performance evaluation. 


A similar situation facing the decision maker is the selection of 
alternative designs of hospitals, schools, housing complexes, or research 


laboratories. 
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Another area of interest is compensation programs to provide employees 
many different options. Companies who can provide many options will 


have a definite advantage in employee inducement and retention. 


4. Long-range large scale programs comprised of many projects that are 


interrelated and interdependent. 


Many government programs such as health care, welfare, education, and 
foreign aid programs are in this category. This is an area of great 
importance because of its far-reaching effects, both intended and 
unintended. It is also the area of greatest difficulty because of its 
complexity, unclear boundary (sphere of effects are not clearly definable) 


and future-oriented consideration. 


These programs or measures that tend to create many side effects or that 
produce long-term effects or irreversible conditions, must receive extra 
care in planning. Although the future is always uncertain and, therefore, nc 
forecasting techniques can claim total accuracy, a variety of forecasting 
techniques combined with modeling can produce some indication of types of 


impacts a given program might make in the future. 


After possible consequences of alternative courses of action are generated, 
the consequences can be arranged within a "decision-event map," indicating 
interrelation of actions taken, their intended results and possible side 


effects, and intervening events that are likely to happen. 
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Concentrating on the consequences in the time-stream (rather than at 

one point in time), our technique can still be employed by using attri- 
butes that explicitly indicate future impacts. (e.g., "rate of yearly 
increase in food production in country X, during 1970-1975, after 
introduction of farming equipment" or "number of farm workers in country 


X migrating yearly into cities during 1970-1975"). 


Admittedly, any future-oriented evaluation is very tenuous. However, 
evaluating programs to assist underdeveloped countries is a more 
amenable problem than evaluating our own future possibilities. We can 
use the U.S. and other developed countries as models in planning to 
avoid possible undesirable consequences and to promote those attributes 
that are desired by the country. Although exact correspondence between 
the model and the real consequences in a given country cannot be 


expected, hindsight is readily available while foresight is not. 
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SUMMARY _ 

The importance of including many criteria of various types and degrees 
of imprecision has been discussed. The man-machine fuzzy-set approach 
described here is our first attempt to tackle this task. Insight gained 
in using these techniques may lead to improvements or to new ideas and 


techniques. 


Systematic analyses of the situation supplemented by intuitive judgment was 


emphasized. The following points may be worth reviewing: 


‘ Consistency in treatment of all alternatives with many attributes 
describing desirability or undesirability. One aspect of consis- 
tency achieved here is the making of everything into a value- 
oriented judgment; even though attribute values may be factual, 
determining their worth in relations to the objectives requires a 
judgmental decision. The fuzzy-set concept allows explicit 


treatment of imprecise value judgments. 


; Comparability. Since absolute judgment is far more difficult than 
relative judgment, the man-machine techniques facilitate comparison 
by bringing in other relevant factors. In addition, the fuzzy-set | 
treatment of attribute values make them commensurable, and complex 
trade-off possibilities can be explored much more readily than 


without such assistance. 
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Systematic use of the knowledge and experience of experts as well as 
opinions of people from different backgrounds. Those techniques 
(such as the Delphi) for direct involvement of people can fit 


naturally into the on-line interactive system. 
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